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SUMMARY 


This  report  describes  two  computer  programs  devel¬ 
oped  to  analyze  the  response  of  multi-cell  rectangular 
box-structures  to  external  dynamic  loads.  The  first  of 
the  two  programs  considers  the  dynamic  response  of  the 
box-structure  with  its  material  properties  remaining 
linear.  The  output  consists  of  displacements,  moments 
and  stresses.  Also,  response  spectra  may  be  generated 
at  equipment  attachment  points  within  the  structure. 

The  second  program  calculates  the  dynamic  response  with 
the  added  provision  that  when  moments  in  any  plate  ex¬ 
ceed  a  predetermined  yield  value,  that  plate  will  not 
accept  any  increase  in  loads  normal  to  its  plane. 

The  report  applies  to  shallow  buried,  reinforced 
concrete  protective  structures.  The  walls,  floors, 
roof  and  partitions  of  these  structures  are  assumed  to 
be  mutually  connected  at  right  angles.  The  structure 
can  have  more  than  one  room  in  the  width,  depth  and 
height.  The  loading  is  that  caused  by  an  above  ground 
nuclear  detonation  some  distance  from  the  location  of 
the  structure. 

The  method  of  analyzing  the  structure  is  by  con¬ 
sidering  it  as  a  system  of  finite  elements,  each  plate 
being  divided  into  four  rectangular  elements.  The 
finite  element  net-hod  is  a  means  of  approximating  the 
equations  of  equilibrium  of  the  elastic  continuum  with 
a  set  of  algebraic  ~quilibriura  equations.  The  continu¬ 
um  is  considered  as  being  an  assemblage  of  discrete 
structural  elements  inter-connected  at  a  finite  number 
of  nodal  points.  The  analysis  involves  the  evaluation 
of  the  element  elastic  properties  which  are  represented 
by  the  stiffness  matrix  expressing  the  relationship 
between  element  nodal  forces  and  displacements*  The 
nodal  points  are  taken  as  the  four  comers  of  the  ele¬ 
ment. 


The  basic  operation  in  the  definition  of  an  element 
stiffness  matrix  is  the  choice  of  deformation  character¬ 
istics  which  are  to  be  allowed.  The  most  important 
criterion  to  be  considered  is  the  compatabilit)  of  de¬ 
formations  of  adjacent  elements.  The  element  stiffness 
matrix  has  twenty-four  degrees  of  freedom,  six  per  nodal 
point,  and  consists  of  a  12x12  bending  stiffness  matrix 


and  a  12x12  inplane  stiffness  matrix.  The  bending 
stiffness  is  calculated  from  a  twelve  term  polynomial 
for  the  displacement.  The  inplane  stiffness  consists 
of  an  8x8  linear  force-displacement  matrix  and  a  4x4 
moment-rotation  matrix.  For  the  purpose  of  this  work 
these  were  assumed  to  be  uncoupled.  At  any  nodal 
point  it  is  possible  to  have  elements  intersecting  in 
three  mutually  perpendicular  planes.  This  makes  it 
necessary  to  convert  the  local  coordinates  for  each 
element  into  a  global  coordinate  system  for  the  whole 
structure. 

At  any  instant  of  time,  the  applied  distributed 
load  on  the  structure  can  be  specified.  This  is  con¬ 
verted  into  a  nodal  point  loading  equivalent  to  the  dis¬ 
tributed  load  such  that  the  nodal  deformations  and  the 
total  work  done  are  the  same  in  each  system.  By  using 
a  process  similar  to  that  used  in  determining  the  stiff¬ 
ness  matrix  for  the  complete  structure,  the  load  matrix 
for  the  complete  structure  can  be  obtained  by  taking  the 
load  matrices  for  each  individual  plate  and  expanding 
them  into  a  global  coordinate  system. 

Through  the  finite  element  method,  the  structure  is 
replaced  by  lumped  parameters  at  discrete  nodal  points, 
therefore,  the  mass  which  is  distributed  throughout  the 
structure  is  assumed  to  be  concentrated  at  those  nodal 
points.  In  order  to  get  the  response  of  the  structure 
for  a  time-dependent  load,  the  equilibrium  of  the  dis¬ 
crete  mass  system  at  time  t  is  expressed  by  the  matrix 
equation 

mro*  ♦  tcux)t  ♦  -  {P)t 

where  CM]  ■  mass  matrix 

[C]  «  damping  matrix 

[K]  -  stiffness  matrix 

(x)t  *  displacement  of  the  system 

{P}^  »  force  acting  on  the  system 


The  equations  are  solved  by  means  of  a  3tep  by  step 
procedure.  The  accelerations  of  the  masses  are  assumed 
to  be  linear  in  each  time  increment.  A  diagonal  mass 
matrix  is  used.  Since  the  determination  of  mode  shapes 
and  frequencies  is  not  a  part  of  the  step  by  step  method 
of  solution,  modal  damping  cannot  be  used  directly.  The 
damping  matrix  is  assumed  to  be  a  linear  function  of  the 
mass  matrix  and  the  stiffness  matrix.  The  solution  of 
the  equations  makes  use  of  the  Gaussian  elimination  tech¬ 
nique. 

An  example  of  a  dynamic  analysis  is  given. 

A  synoptic  review  of  selected  literature  pertaining 
to  the  failure  modes  and  failure  loadings  of  two-way 
slabs  is  presented  at  the  end  of  the  report. 


PREFACE 


During  the  year  of  1963*  the  Office  of  Civil  De¬ 
fense,  Protective  Structures  Division,  initiated  the 
development  of  a  series  of  computer  programs  directed 
toward  providing  a  complete  engineering  analysis  of 
buildings  considerin':  the  multiple  effects  of  nuclear 
weapons.  These  computer  codes  are  the  backbone  of  the 
Computer  Analysis  For  Protective  Structures  (CAPS) 
system. 

The  CAPS-1  program,  entitled  Analysis  of  Structures 
for  Resistance  to  Nuclear  Blast ,  began  with  the  devel¬ 
opment  of  a  code  for  determining  the  blast  resistance 
of  multi-3tory  buildings.  T.Y,  Lin  and  Associates  com- 

Sleted  the  computer  program  under  Contract  No.OCD-OS-63- 
4,  Subsequently,  the  programming  effort  to  provide  a 
code  for  analysing  box-structures  was  begun  under  Con¬ 
tract  No.  OCD-PS-&5-7,  Subtask  1157A.  The  responsibility 
for  continuation  of  the  program  passed  to  Stanford  Re¬ 
search  Institute  under  the  provisions  of  Contract  No, 
OCD-PS-6tl-201  between  OCD  and  SRI.  T.  Y.  Lin  and 
Associates  has  performed  the  work  contained  in  this  re¬ 
port  under  Subcontract  No.B-8l869(*l9,*9A-A6)-US  with  SRI, 
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INTRODUCTION 
1-1  Background 

This  report  describes  the  initial  development  of  a 
computer  program  designed  to  analyze  the  dynamic  response 
of  box-type  structures  to  nuclear  blast.  A  box-type 
structure  generally  carries  external  loads  by  3labs 
spanning  between  bearing  walls  or  shear  walls.  When 
used  as  shelters,  these  structures  may  be  aboveground, 
located  in  the  basement  of  a  larger  building,  or  burled, 

A  previous  report,  by  T.  Y.  Lin  and  Associates1* 
under  Contract  No.  OCD-PS-65-7*  drew  together  pertinent 
conclusions  of  a  large  number  of  reasearch  papers  and 
design  manuals  in  order  to  bring  them  to  bear  on  the 
development  of  this  computer  program.  Airblast  loadings 
and  ground  motions  were  reviewed.  Also,  the  equations 
for  the  dynamic  response  of  lumped  mass  systems  were 
presented  in  detail. 

Many  other  publications  have  presented  design- 
analysis  methods  applicable  to  boxes2»3>%»5.  In  general, 
an  element-by-eleraent  approach  is  used.  That  Is,  the 
loads  are  applied  to  the  primary  resisting  element,  the 
response  is  determined  from  a  single-degree-of-freedom 
idealization,  and  the  reactions  subsequently  applied  as 
loads  to  the  supporting  members. 

For  purposes  of  this  computer  program,  the  box  has 
been  thought  of  as  a  multl-degree-of-freedom  system 
consisting  of  plates,  rigidly  connected  to  each  other 
at  right  angles  along  each  edge.  The  response  analysis 
is  a  finite  element  technique,  utilizing  rectangluar 
plate  elements. 

Other  publications  have  described  al:  Mast  loads  on 
an  aboveground  box 2,3 •  Loads  on  an  underground  box  are 
not  well  defined.  The  fundamental  phenomena  of  the  soil- 
structure  interaction  process  are  not  completely  under¬ 
stood.  However,  the  relative  compressibility  of  the  soil 
and  structure  probablv  controls  the  load  experienced  by 
the  buried  structure, 


•Superscripts  refer  to  references  found  in  Section  7. 
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The  computer  program  models  the  compressibility  of 
a  box.  An  adequate  model  of  the  soil  is  still  needed  be¬ 
fore  an  analysis  of  a  soil-structure  system  can  be  under¬ 
taken  with  any  degree  of  confidence. 

1-2  Scope  of  the  Report 


The  program  development  is  divided  into  four  major 
phases. 

Phase  I  Development  of  a  static  analysis  for 

box-type  structures  utilizing  rectangu¬ 
lar  finite  elements 


Phase  II  Extension  of  the  static  analysis  to 

elastic  dynamic  response  analysis  using 
a  step-by-step  integration  procedure 

Phase  III  Incorporate  a  procedure  for  determining 
the  response  spectrum  at  each  nodal 
point  of  the  elastic  system 

Phase  IV  Incorporate  non-linear  definition  of 
plate  bending  property  into  Phase  II 


This  report  covers  the  complete  mathematical  formu¬ 
lation  of  each  phase.  An  example  of  an  elastic  multi¬ 
cell  box  analysis  for  dynamic  loads  is  given.  An  example 
of  a  non-linear  single  cell  box  analysis  for  dynamic 
loads  is  also  given,  in  which  plates  fail  successively 
as  they  exceed  a  given  yield  capacity. 


The  dynamic  loads  are  assumed  to  act  on  the  exterior 
surfaces  of  the  multi-cell  box.  Time  did  not  permit  the 
programing  of  a  comprehensive  load  routine.  However,  a 
simple  formulation  of  loads  on  a  shallow  buried  box  was 
added  to  the  main  program  in  order  to  carry  out  an  example 
analysis  of  a  box  responding  to  a  time-dependent  load. 

The  load  subroutine  is  described  in  Appendix  A. 
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DEFINITION  OF  THE  STRUCTURE 
2-1  General 

Many  underground  protective  structures  are  built  of 
heavy  concrete  walls,  with  a  concrete  floor  and  a  thick 
concrete  roof.  The  structures  may  have  more  than  one 
level  in  which  case  the  intermediate  floors  are  concrete 
also.  Thus,  the  whole  structure  can  be  considered  as  a 
set  of  interconnecting  plates  mutually  at  right  angles; 
l.e.,  a  box-type  structure* 

The  analysis  of  plates  at  right  angles  presents  two 
major  problems.  The  first  is  that  two  systems  of  coord¬ 
inates  have  to  be  used,  l.e.,  a  local  system  for  the  in¬ 
dividual  plates  and  a  global  system  for  the  whole  struc¬ 
ture.  Extreme  care  must  be  taken  when  changing  from  one 
to  the  other. 

The  second  major  problem  is  the  large  number  of 
mutually  dependent  coordinates.  In  the  case  of  three 
dimensional  frame  analyses  the  numbering  of  dependent 
coordinates  may  be  arranged  in  such  a  manner  that  the 
stiffness  matrix  can  be  formulated  with  a  narrow  band 
width;  thus,  core  space  in  the  computer  is  conserved. 

In  the  case  of  plates  mutually  connecting  at  right  an¬ 
gles,  the  band  width  can  become  very  large,  thereby 
creating  a  storage  problem. 

An  accurate  method  of  analyzing  box-type  structures 
and  the  method  adopted  in  this  work,  is  to  define  the 
box  as  a  system  of  finite  elements.  Each  plate  can  be 
considered  as  an  element  or  further  subdivided  Into  a 
number  of  elements.  The  essential  feature  of  the  finite 
element  method  is  the  means  by  which  the  differential 
equations  of  equilibrium  of  the  elastic  continuum  are 
approximated  by  a  set  of  algebraic  equilibrium  equations. 
This  procedure  is  generally  looked  upon  as  the  substitu¬ 
tion  for  the  actual  continuum  of  an  assemblage  of  dis¬ 
crete  structural  elements.  Interconnected  at  a  finite 
number  of  nodal  points.  In  effect,  the  continuum  may 
be  visualized  as  being  physically  cut  up  into  the  finite 
element  system,  the  material  properties  of  the  original 
material  being  retained  in  the  elements.  The  analysis 
involves  the  evaluation  of  the  element  elastic  prop¬ 
erties,  which  are  represented  by  the  stiffness  matrix 
expressing  the  relationship  between  element  nodal  forces 
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and  displacements.  The  nodal  points  are  taken  as  the 
four  corners  of  the  element. 

With  this  system,  it  is  possible  to  define  each 
panel  as  an  element,  determine  the  equivalent  nodal 
point  stiffness  and  load  and  solve  for  displacements. 

From  the  displacements,  the  moments  and  shears  at  the 
nodal  points  can  be  determined.  In  order  to  determine 
the  distribution  of  the  moments  and  shears  throughout 
the  panels,  it  is  necessary  to  subdivide  each  panel  in¬ 
to  smaller  elements.  The  smallest  number  of  subdivisions 
to  obtain  the  moments  and  shears  at  the  center  of  each 
panel  is  four  elements.  As  the  least  complicated  box, 
l.e.,  a  cube,  has  six  sides,  it  is  apparent  that  the 
stiffness  matrix  can  soon  get  impossibly  large.  For 
this  work  it  was  decided  that  a  four-element  panel  was 
the  least  that  could  be  used  to  get  results  for  moments 
and  shears  in  the  middle  of  the  panels  and  that  larger 
subdivisions  would  be  unworkable. 

It  was  realized  that  the  development  of  the  stiff¬ 
ness  matrix  for  the  structure  in  the  global  coordinate 
system  would  require  some  difficult  programming,  so  the 
preliminary  work  was  done  on  a  cube.  A  symmetrical 
static  load  wa3  applied  about  one  major  plane  so  that 
symmetrical  displacements  would  constitute  one  check  on 
the  accuracy  of  the  stiffness  matrix. 

When  the  stiffness  matrix  was  considered  correct, 
steps  were  taken  to  convert  the  program  to  a  system  of 
plates  divided  into  four  elements.  The  correctness  of 
this  change  was  determined  by  comparing  the  analysis  of 
a  one-box  structure  using  the  four-element-per-panel 
system  with  the  same  structure  modelled  as  an  eight-box 
structure  using  the  one-element-per-panel  system.  In 
the  latter  case,  the  stiffness  of  the  interior  elements 
was  set  to  zero. 

The  method  of  determining  the  nodal  point  loading 
for  a  time  dependent  blast  force  was  solved  In  the  early 
stages  of  the  work.  By  considering  this  load  at  any  given 
time  as  a  statically  applied  load,  it  was  possible  to  get 
representative  moments  and  shears  in  the  structure.  The  , 
procedure  adopted  resulted  in  a  method  of  complete  check¬ 
ing  as  every  stage  was  completed  and  provided  for  a 
logical  flow  of  work  from  a  system  of  single-element 
panels  under  static  loads  to  a  system  of  four-element 
panels  under  dynamic  loads. 
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2-2  The  Stiffness  Matrix 

The  box-type  structure  is  considered  as  a  set  of 
interconnecting  plates  in  three  mutually  perpendicular 
planes.  In  order  to  determine  the  stiffness  matrix  for 
the  complete  box,  it  is  necessary  to  determine  the 
stiffness  of  each  plate  in  its  own  local  coordinate 
system  and  then  convert  this  to  the  global  system.  The 
definition  of  these  systems  is  shown  in  Figure  2-1. 

The  development  of  the  system  for  considering  four 
elements  per  panel  is  essentially  the  same  as  that  for 
the  system  with  one  element  per  panel  with  the  one  dif¬ 
ference  that  there  will  be  restraining  panels  at  right 
angles  on  only  two  of  the  edges.  For  simplicity  of  pre¬ 
sentation,  the  stiffness  matrix  of  the  structure  will  be 
developed  as  if  each  panel  consists  of  only  one  element. 

Since  the  method  used  is  to  consider  each  plate  as 
a  finite  element  in  the  box,  by  assuming  Klrchoff  plate 
bending  theory,  the  stiffness  is  determined  in  terms  of 
the  24-comer  displacements  of  each  plate.  The  inter¬ 
section  of  the  comers  of  the  plates  are  the  nodal  points 
of  the  complete  system. 

The  basic  operation  in  the  definition  of  an  element 
stiffness  matrix  is  the  choice  of  the  deformation  charac¬ 
teristics  which  are  to  be  allowed.  The  most  important 
criterion  to  be  considered  in  this  selection  is  the  com- 
patabllity  of  deformations  of  adjacent  elements.  It  can 
be  shown  that  if  deformation  prcteros  are  specified  which 
provide  internal  compatability  within  the  elements,  and 
at  the  same  time  achieve  full  compatability  of  displace¬ 
ments  along  the  element  boundaries,  then  the  strain 
energy  in  the  idealization  will  represent  a  lower  bound 
to  the  strain  energy  of  the  actual  continuum. 

The  relationship  between  the  externally  applied 
forces  and  the  resulting  displacement  on  a  plate  in  its 
own  local  coordinate  system  can  be  expressed  in  the 
form 


(P>  ■ 


(u) 
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where  S,  is  the  12x12  bending  stiffness  matrix 
3-  is  the  12x12  inplane  stiffness  matrix. 

There  is  t*Io  coupling  between  these  two  matrices. 

The  simplest  expression  which  has  been  used  in  de¬ 
fining  the  rectangular  element  bending  stiffness  is 
the  12-term  nolynominal  in  x  and  t/;13 

10  a  a3  ♦  a2X  +  a3t/  ♦  «i,x2  +  as*(t  +  a6 ^  +  a7*3 

*  «S2y  *  a«xy2  ♦  «iQy3  +  an*V  +  ai2*y3  (2_v 


The  two  fourth  degree  terms  which  are  included  to 
provide  the  required  number  of  displacement  patterns  were 
selected  to  maintain  symmetry  and  also  to  satisfy  the  bi¬ 
harmonic  equilibrium  equation.  A  complete  12  by  12 
stiffness  matrix  which  represents  the  nodal  force  deflec¬ 
tion  relationships  is  shown  in  Figure  2-2,  This  12  by  12 
matrix  may  be  constructed  from  the  12  by  3  matrices  asso¬ 
ciated  with  each  Joint.  For  any  given  joint  x,  the 
stiffness  matrix  can  be  represented  as  the  sum  of  four  12 
by  3  matrices  KA,  *b»  Kc  and  Kd ,  i.e.,  •  Ka  *  Kb  +  *c 

♦  Kd*  These  represent  the  energy  contributions  from  the 
("xx)2.  (wi/</)2»  (vxx*wytf)  and  (*xu)2  curvature  terms,  re¬ 
spectively.  The  coefficients  of  these  four  stiffness 
component  matrices  as  developed  by  Adini7  are  shown  in 
Figure  2-3. 


The  12  x  12  inplane  stiffness  matrix  can  be  con¬ 
sidered  as  an  8  x  8  linear  force-displacement  matrix  S. 
and  a  4  x  moment-rotation  matrix  Sr  together  with  the 
corresponding  coupling  terms  and 


St 

- -  I—  -  * 

l  ri\r  . 


In  order  to  determine  the  S^r  and  Sp£  terms  a  rigorous 
solution  la  required,  according  to  the  theory  of  elasticity. 
The  authors  have  assumed  the  accuracy  requirements  do  not 
warrant  this,  and  expressions  for  the  linear  and  rotational 
displacements  cere  developed  as  if  uncoupled. 

For  the  8x8  linear  stiffness  matrix  S.  the  dis¬ 
placement  pattern  is  as  shown  in  Figure  2-A.  The  sides 
remain  straight  after  deformation. 
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Note 
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Pi sure  2-2  ARRANGEMENT  OP  BENDING  STIPPNESS  MATRIX 
(Reference  13  -  Note:  All  minus  signs  omitted  in  reference) 


Figure  2-3  ADIMI  COEFFICIENTS  (Reference  7) 
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-{uy  -  U£> 
£l“k  "  U4) 
HF(u-t  '  tty  ' 


or  in  matrix  form, 


I  -I  -1 


(214) 


the^same^result^l  ,e!'°Ce^Ure  f°r  th*  “  diSpla0er,e"ts  »!•*»• 


b  0  0 

b  b  0 
a  0  a 
I  -1  -1 


(2-5) 


From  the  definition  of  strain; 

«n 


Hence,  the  relationship  between  the  displacements  and  the 
strains  can  be  expressed  in  the  form 

1  c x  [-fa*,,  6-y  -y  y  0  0  0  0  ' 

;*»  mh  *  *  0  » -«♦«  -*  «-«  *  !  «,  j 

y  !  -a*x  -x  a-x  x  -£>*</  b-t?  -yy  j 
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The  rotational  stiffness  matrix  is  difficult  to  de¬ 
fine.  Because  of  this  fact,  the  procedure  in  plate  and 
shell  theory  has  been  to  ignore  it.  For  the  box  type 
structure  this  is  not  possible  so  displacement  functions 
were  decided  on  which  would  represent  the  rotation  of 
the  corners  in  turn. 


y 


INPLANE  ROTATIONAL  STRAIN  OF  AN  ELEMENT 


The  displacement  pattern  is  shown  in  Figure  2-5  where 
u(x,y)  -  (</-i>)2U-x)  -  (2-7) 

v(x,</)  *  U-<02(fc-</)  *  F  (x,</)e  (2-8) 

a2b  2 
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Since  the  assumed  displacement  functions  have  a 
linear  0  dependence,  superposition  can  be  used. 

Hence 

u(x ,{f)  *  “u)  *  tt(yj  *  u(fe}  *  uU) 
v{K,y)  •  vu)  ♦  w(yj  *  v(fe)  ♦  ^u) 

or  in  matrix  form  as  shown  in  figure  2-6. 

The  relationship  between  the  displacements  and  the 
strains  can  be  exnressed  in  the  form 


or 

[e]  =  [a  ]{e}  (2-9) 

0 

where  [a0]  is  given  in  Figure  2-7. 

Due  to  the  displacement  field  (u,v)  for  the  set 
of  angle  changes  (e)  there  Is  a  field  of  external  self 
equilibrating  moments  [M], 

The  equilibrium  condition  is  satisfied  if  the  set 
of  external  moments,  [M],  acting  on  the  external  angle 
changes,  {8}  produces  the  same  work  as  the  set  of  in¬ 
ternal  stresses,  {o},  acting  on  the  internal  strains, 

u>. 

Hence  Me>T[M]  -  /  U  >T  (o  }d{vof )  (2-10) 

vol 

Using  equations  2-6  ana  2-9 

ab 

(0>TCM]  -  <J/(e}r[«eJr[K3[aJ(e>c(xd</  (2-U) 

oo  B 

Since  (8)  is  independent  of  "x*  and  mifm  , 

CM]  -  [c//ta85TCK3Cae]c(x<iS{Jt8j  <2-l« 

00 
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Figure  2-6  ASSUMED  DISPLACEMENT  FUNCTION 
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Figure  2-T  DISPLACEMENT  STRAIN  RELATIONSHIP 
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So  that, 


[SI  -  ^//C<t„]T[K]C<t  Idxdg  <2-13) 

00  9  0 

The  expansion  of  the  integral  in  equation  2-13  for 
the  first  row  of  the  matrix  Cs_ ]  is  shown  in  Figure  2-8, 

It  can  be  seen  that  there  is  a  lack  of  equilibrium, 
therefore  the  matrix  was  brought  into  equilibrium  by 
arbitrarily  subtracting  from  each  element.  This  is 
numerically  small  compared  with  the  other  terms  and 
ha3  little  effect  on  the  results. 

The  first  row  of  the  matrix  C 3  is  expanded  into 
a  consistent  matrix  as  shown  in  FigSre  2-9. 

As  plates  in  different  planes  are  joined  together, 
the  elements  of  the  stiffness  matrices  cannot  be  added 
directly.  It  is  necessary  to  convert  the  coordinates 
from  a  local  to  a  global  system.  This  can  be  done  by 
using  the  coordinate  transformation  matrix  for  stiffness 
as  shown  in  Table  2-1 .  Once  in  the  global  system  the 
stiffness  matrices  can  be  combined  into  a  single  stiff¬ 
ness  matrix  [K]  for  the  complete  box  by  keeping  track  of 
the  nodal  points  of  the  four  corners  of  the  individual 
elements. 
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TABLE  2-1 


Coordinate  Transformation 
for 

Stiffness  Coefficients 


LOCAL 

I 

GLOBAL 

I  ^ 

Type  1  Type  2 

Type  3 

1 

1 

3 

3 

2 

7 

9 

9 

3 

13 

15 

15 

4 

19 

21 

21 

5 

-2 

-2 

-1 

6 

-8 

-8 

-7 

7 

-14 

-14 

-13 

8 

-20 

-20 

-19 

9 

-3 

1 

-2 

10 

5 

5 

4 

11 

4 

6 

6 

12 

-9 

7 

-8 

13 

11 

11 

10 

14 

10 

12 

12 

15 

-15 

13 

-14 

16 

17 

17 

16 

17 

16 

18 

18 

18 

-21 

19 

-20 

19 

23 

23 

22 

20 

22 

24 

24 

21 

-6 

4 

-5 

22 

-12 

10 

-11 

23 

-18 

16 

-17 

24 

22 

-23 
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Section  3 


MODAL  POINT  LOADS 

At  any  instant  of  time,  the  applied  distributed  load 
on  the  structure  can  be  specified.  This  must  be  converted 
into  a  nodal  point  loading  equivalent  to  the  distributed 
load  such  that  the  nodal  deformations  and  the  total  work 
are  the  same  in  each  system. 

h  w[xty)dA  (3'1.) 

A 

hU)T  U>  (3-2) 


where  q  is  the  load  on  the  area  <M  and  to  is  its  displace¬ 
ment.  For  one  plate,  U)  is  a  12  x  1  column  matrix  con¬ 
sisting  of  a  force  and  two  moments  at  each  corner  of  the 
plate,  and  {$}  is  the  associated  displacement. 

From  Equations  3-1  and  3-2, 


Substituting  the  12-term  polynomial  for  to  in  Equation  2-1, 

U)rU)  -  |  U3U}tM 

where  [tf]  is  a  1  x  12  matrix  of  the  terms  of  the  polynomial 
and  (i)  Is  a  12  x  1  matrix  of  the  coefficients. 

Substituting  for  [a]  • 

U,T-  l 


-  CN]r  |  rrf* 


or 

(U 


External  Energy 
due  to  actual  * 
distributed  load 

External  Energy 
due  to  equivalent  * 
nodal  point  loads 
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The  matrix  [M]  which  has  been  introduced  can  be 
evaluated  through  the  polynomial  in  w  with  this  expres¬ 
sion  and  by  differentiating  it  with  respect  to  x  and  yt 
the  values  of  10  can  be  related  for  each  displacement  at 
each  corner  for  the  twelve  coefficients  o.  The  develop¬ 
ment  of  the  inverse  matrix  Cm3”1  is  given  in  Figure  3-1. 

The  function  can  be  expressed  as: 


«(x#0)  *  C*  +  C2X  +  C3 !/  +  CHxy 


where  the  coefficients  CJ#  C2,  C,,  and  C.  are  evaluated 
from  the  load  conditions  at  the  nodal  points. 

The  integral  /  q.  .[d]^rfA  can  be  written  as: 


U)  =*  f(Cl  +  C2r  +  C  3  (/  +  C4  xi/)  [d]^dA 
A 


or 

(I)  =  [F(tt,b 

r  'i 

„  -1 

and 

[$] 

4 

qj 

►  _  4 

«k 

\ 

qt 

V 

Cl 

c 

C3 

C„ 


so  that  U)  »  [M]rCF(a#6|  ][!?]  {<?} 


(3-5) 


where 


C  <T‘] 


100  0 

1  a  0  0 

10b  0 

I  a  b  ab 


and 


is  shown  in  Figure  3-2. 
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Figure  3-1  DKVT.1,01  .'I-'>JT  OF  MATRIX  [N ] 
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MATRIX  F(a,b) 


By  using  a  procedure  similar  to  that  used  in  deter¬ 
mining  the  stiffness  matrix  for  the  complete  structure, 
the  load  matrix  for  the  complete  structure  can  be  obtain¬ 
ed  by  taking  the  load  matrices  for  each  individual  plate 
and  expanding  them  into  a  global  coordinate  system  by  us¬ 
ing  a  coordinate  transformation  matrix  and  inserting  the 
elements  of  the  submatrices  in  the  large  matrix  at  the 
desired  location  .  The  load  coordinate  transformation 
matrix  is  shown  in  Table  3-1. 


TABLE  3-1 


Coordinate  Transformation  for  Nodal  Loads 


LOCAL 

GLOBAL 

Type  1 

Type  2 

Type  3 

1 

-3 

1 

-2 

2 

5 

5 

4 

3 

4 

6 

6 

4 

-9 

7 

-8 

5 

11 

11 

10 

6 

10 

12 

12 

7 

-15 

13 

-14 

8 

17 

17 

16 

9 

16 

18 

18 

10 

-21 

19 

-20 

11 

23 

23 

22 

12 

22 

24 

24 
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Section  *1 

DYNAMIC  RESPONSE  ANALYSIS  -  ELASTIC 
*1-1  Equation  of  Equilibrium 

Through  the  finite  element  method,  the  box  structure 
is  replaced  by  lumped  parameters  at  discrete  nodal  points. 
The  mass,  which  is  distributed  throughout  the  structure, 
is  assumed  to  be  concentrated  at  the  nodal  points.  The 
deformation-resistance  characteristics  of  the  box  are  de¬ 
fined  for  relative  displacement  of  the  mass  coordinates  by 
the  stiffness  matrix  described  in  Section  2, 

The  equilibrium  of  the  discrete  mass  system  at  time 
t  is  expressed  by  the  following  matrix  equation: 

CM]  {X>^  +  [C]{X>^  +  [K]  {X)^  =»  (P)t  (4-1) 

•  * 

where  {X}^  »  acceleration  of  the  system 
(XI  *  velocity  of  the  system 
{X}^  »  displacement  of  the  system 
{P}^  »  force  acting  on  the  system 
CM]  »  the  diagonal  mass  matrix 
[C]  *  the  damning  matrix 
(X]  ■  the  stiffness  matrix 
4-2  Mass  Matrix 

A  lumped  mas3  approach  la  used  to  calculate  the  dynam¬ 
ic  response,  with  the  displacement  of  each  mass  directly 
represented  by  a  single  coordinate.  The  advantage  of  this 
simple  mass  idealization  is  that  it  eliminates  any  mass 
coupling.  Thus,  a  diagonal  mass  matrix  is  used  in  enuation 
4-1. 


The  mass  at  each  nodal  point  of  the  structure  con¬ 
tains  a  contribution  from  each  plate  element  attached  to 
the  node.  Referring  to  Figure  4-1,  thi3  typical  mass 
contribution  at  a  node  from  a  single  element  is  expressed 
as, 
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where  p  and  t  are,  respectively,  the  density  and  thick¬ 
ness  of  the  plate  element.  Similarly,  p  and  t8  are, 
respectively,  the  density  and  thickness  of  soil  layer 
assumed  acting  with  the  plate  element. 

In  equation  4-2,  the  distributed  mass  of  each  plate 
element  i3  assumed  lumped  at  its  corners,  one-fourth  of 
the  total  mass  to  each.  Figure  4-1  illustrates  this 
physical  lumping  of  mass.  The  mass  at  each  corner,  or 
nodal  point,  can  be  put  into  motion  in  the  local  x,  y  and 
z  directions.  Mass  quantities  associated  with  these  move¬ 
ments  are  expressed  by  the  first  term  of  the  right-hand 
side  of  the  equation.  This  mass  matrix  In  its  diagonal¬ 
ized  form  is  an  extension  of  the  w ork  by  Archer  on  con¬ 
sistent  mass  matrices  for  distributed  mass  systems.12 

Since  the  structure  is  in  contact  with  soil,  many 
exterior  plate  motions  wJll  be  coupled  with  the  movement 
of  a  certain  mass  of  soil.  The  second  term  of  Equation 
4-2  assigns  a  mass  quantity  to  the  local  z  direction  of 
element  displacement  only.  In  other  words,  soil  inertia 
forces  are  considered  only  for  motions  perpendicular  to 
each  element's  plane.  The  thickness  of  soil,  is  an 
arbitrary  quantity.  For  the  roof  of  a  shallow  buried 
box,  f  equals  the  depth  of  soil  cover. 
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Figure  *4-1 


MASS  CONCENTRATION  AT  A  TYPICAL  NODE 


TABLE  *1-1 


GLOBAL 

LOCAL 

Type  1  Type  2  Type  3 

1 

1  3  3 

2 

2  2  1 

3 

3  1  2 

*i 

k  6  6 

5 

5  5  * 

6 

6  *i  5 

COORDINATE  TRANSFORMATION  FOR  MASSES 


*»-3 


The  mass  matrix  for  the  whole  structure  is  obtained 
by  expanding  the  mass  matrix  for  each  individual  element 
into  the  global  coordinate  system  by  means  of  a  coordi¬ 
nate  transformation,  and  superposing  the  expansions. 

A  generalized  mass  approach,  based  on  assumed  dis¬ 
placement  functions,  could  be  used.  However,  this  would 
lead  to  coupled  (i.e.,  off-diagonal)  mass  matrix  terms 
which  complicate  the  storage  problems  in  the  computer. 
Furthermore,  assumed  shape  functions  will  not  necessari¬ 
ly  lead  to  improved  accuracy  in  the  results.  The  most 
severe  approximation  is  the  number  of  degrees  of  freedom 
allowed.  It  is  believed  that  the  physical  lumping  of 
mass  is  a  satisfactory  approximation  for  this  particular 
system. 


4-3  Damping  Matrix 

The  exact  form  of  the  damping  matrix  is  unknown  for 
box-structures,  as  is  the  case  for  most  structural  systems. 
However,  in  most  instances,  its  effect  on  the  mode  of  vi¬ 
bration  will  be  small.  Therefore,  it  i3  reasonable  to 
replace  the  [C]  matrix  by  the  following  matrix  relation¬ 
ship:  1 


CC]  -  ao[M]  +  ^[K]  (4-3) 

By  assuming  the  damping  matrix  to  be  a  linear  combination 
of  the  mass  matrix  and  the  stiffness  matrix,  computations 
are  simplified  and  computer  storage  requirements  are  mini¬ 
mized. 

The  significant  test  experience  in  determining  energy 
losses  due  to  structural  damping  has  related  these  losses 
to  the  frequencies  and  mode  shapes  of  particular  systems. 
Since  the  determination  of  frequencies  and  mode  shapes  is 
not  a  part  of  the  step-by-step  method  of  solution,  modal 
damping  cannot  be  used  directly.  However,  «0  and  «  may 
be  interpreted  in  terms  of  equivalent  modal  damping1  if 
the  significant  frequency  range  of  response  is  known  or 
can  be  estimated. 

The  modal  damping  ratio  X  •  for  the  Xth  mode  is  given 

in  terms  of  a  and  a  by,  * 
o  | 
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(4-4) 


f  u 

where  Is  the  circular  frequency  of  the  l  mode. 

The  frequency  HT,  which  yields  a  minimum  value  of 
the  damping  ratio  T,  is  given  by 


(4-5) 


If  the  minimum  damping  ratio  X  and  its  associated 
frequency  u  are  given,  the  damping  coefficients  aQ  and  dj 
are  calculated  from  the  following  equations: 


a  *  T  u 
o 

Equation  4-4  can  now  be  rewritten: 


(4-6a) 

(4-6h) 


(4-7) 


Usually,  the  numerical  value  of  the  minimum  damping 
ratio  will  be  established  from  Judgement  regarding  past 
structural  tests.  It  may  vary  through  a  range  of  2  per¬ 
cent  to  30  percent  of  critical  (.<**<  X  <.30).  Having 
estimated  a  minimum  damping,  however,  does  not  fully  de¬ 
fine  the  coefficients  a  and  a^.  The  associated  frequency 
u  oust  be  established.  0 

If  the  significant  frequency  range  of  response  is 
between  the  values  and  t»n,  such  that; 

( 4—8) 

it  is  convenient  to  select  the  frequency  X  as  a  multiple 
of  the  lowest  frequency  and,  in  turn,  such__that  the 
highest  frequency  *  is  the  same  multiple  of  uT. 
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(Ma) 

(4-9b) 


That  I3; 


U  “  n  'll)  j 

w  *  m  w  ■  m2w, 

fl  * 

where  m  i3  any  positive  constant.  Mow,  the  frequency  w 
can  be  calculated  simply  from  the  established  frequency 
range,  as  follows: 


u 


(4-10) 


The  values  of  and  dj  are,  as  given  previously; 

a.  »  T  w  (4-6a) 

0 

-  T/w  (Wb) 


The  damping  ratio  for  the  first  mode,  *  1 ,  and  the 
highest  mode  of  interest,  A  ,  becomes  (rewriting  equation 

4-7); 


(4-11) 


As  an  example  calculation,  assume  the  minimum  damp¬ 
ing  ratio  to  be  5  percent  of  critical.  'Further,  assume 
that  the  fundamental  period  of  vibration  is  associated 
with  the  flexural  response  of  thick  slabs  and  is 
(**i  •  Si.  6  rad/sec.  Also,  assume  that  the  highest  fre¬ 
quency  of  Interest  is  that  due  to  inplane  response  of  the 
slab  elements  and  is  equal  to  316  rad/sec: 
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Z  •  J 31.6  x  3 16  »  100  rad/sec 
T  *  .05 

x  a  ,  m  (31.6  *  316  )  .05  m 
i  * 


a  ■  .05  x  100  "  5.0 
o 

&l  -  .05  ♦  100  -  0.0005 


4-4  Response  Calculations 

It  Is  assumed  that  the  acceleration  of  each  dis¬ 
crete  mass  varies  linearly  within  a  time  interval  A t. 
This  assumption,  which  is  illustrated  in  Figure  4-2, 
leads  to  a  parabolic  variation  of  velocity  and  a  cubic 
variation  of  displacement  within  the  time  interval. 

A  direct  integration  over  the  interval  gives  the 
following  equations  for  acceleration  and  velocity  at  the 
end  of  the  time  interval: 


X  * 

t 

9 

At2 

Xt  *  At 

( 4-l2a) 

*t  • 

5 

At 

Xt  *  8t 

( 4-12b) 

where  A^,  • 

6 

At2 

Xt-At 

6  X  ♦ 

U 

tx 

t-At 

( 4-12c ) 

8*  • 

5 

At 

Xt-At  * 

“t-u  * 

—  Xt-At 
t 

( 4-12d) 
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X* 


Acceleration 


Figure  4-2 

ASSUMED  BEHAVIOR  OF  TYPICAL  DISPLACEMENT  COMPONENT 


4-8 


Substituting  Equations  4-3,  4-12a  and  4-l2b  into  the 
dynamic  equilibrium  relation  4-1  yields  the  following 
set  cf  linear  equations  in  terms  of  some  unknown 
"effective"  displacement. 


03  [X3t 

a 

(4-13) 

where 

03 

a 

[<■]  +  C2tM] 

(4-l4a) 

CP3, 

» 

OV  cjb]! 

(4-l4b) 

3 

I  [X]*  -  Si  [8]. 

Ci  *  c2  t 

(4-l4c) 

in  which 

C 

0 

- 

JL  +•  £^o 

A-t2  u 

(4-15a) 

C, 

M 

1 

(4-l5b) 

i  ♦  3c, 

17” 

C2 

a 

C0C, 

(4-15c) 

- 

C2*l 

(4-15d) 

m 

V  CJ 

(4-15e) 

The  dynamic  response  ^f  the  structure  is  obtained 
from  the  above  relationships  by  using  the  following  step 
by  step  procedure, 

1.  Initialization 

a.  Form  stiffness  matrix  [K]  and  mass  matrix  [/«] 

b.  Form  "effective"  stiffness  matrix 

03  -  03  ♦  C2  O'] 

c.  Trlangularize  03 


4-9 


2,  For  Each  Time  Increment 


a.  Form  [Aj^and  [8]^ 

[A],  -  j-  i-  *[««.** 

Z  A £2  A-t 

[g],  a  JL  ^X^-At+  ?^^-A^+—  tfW 

*  A-£  2 

b.  Form  "effective"  load 

[p]f  »  tp]f  +  Cm3  jm*  +  cjs]t] 

c.  Solve  for  "effective"  displacements 

[x]f  •  [7]-‘  (P]t 

d.  Calculate  Displacements,  Velocities  and 
Accelerations  at  time  £. 


£«t  • 

Ci 

+  Cirt^S] 

£ 

JL 

Af 

[X], 

-  [8], 

•ft 

[XI  * 

6_ 

[X]  . 

-  [A ] 

£ 

u2 

£ 

£ 

s.  Repeat  for  next  time  increment 

It  is  important  to  note  that  for  elastic  structures, 
the  matrix  [X3  need  only  be  formed  once  since  it  is  In¬ 
dependent  of  time. 
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4-5  Solution  of  the  Equilibrium  Equation 

Equation  4-13  expresses  the  dynamic  equilibrium  of 
the  system  in  an  equation  of  the  form,  [A]  [X]  *  [8], 

The  method  used  for  the  solution  of  these  equations 
is  the  Gaussian  Elimination  method. 

4-5.1  Gaussian  Elimination 

The  equilibrium  equations  for  a  system  of  finite 
elements  may  be  written  in  the  following  form: 


X 

♦  A  X 

♦  A  X  .  .  . 

. . .  ♦  A,  X  •  8 

(4-l6a) 

11  1 

12  2 

13  3 

IN  N  i 

X 

21  1 

♦  A  X 

22  2 

♦A  X  ... 
23  3 

...♦  a2Mxm  -  b2 

(4-l6b) 

X 

31  1 

♦  A  X 

32  2 

♦A  X  ... 

33  3 

•••*  A3MXW  *  8 3 

(4-l6c) 

AW1X1  *  A?|2X2  *  AN3X3 


*  ANNXN  “  8N 


The  first  step  in  the  solution  of  the  above  set  of 
equations  is  to  solve  equation  4-l6a  for  Xp  or 


\  '  V\,- 


If  equation  4-17  is  substituted  into  equations  (4-l6a, 
4-l6b,  U)  a  modified  set  of  W-l  equations  is  de¬ 

termined. 

•  A»x» . *  a!mxn  • »!  <*-»«> 
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a*  x  ♦  ALX, 

32  2  33  3 


,♦  A1  X  -  8l 
3M  M  2 


*  *i.X. 


♦  A1  X  »  81 
WW  N  M 


(4-l8b) 


where 


‘i/  •  VA*»VA» 


f  j  *  •  * « •  |  N 


(4-19a) 


•  8^-A^l8l/An  x  -  f, _ _  M 


(4-19b) 


A  similar  procedure  is  used  to  eliminate  X2  from  equation 
4-l8a,  etc. 

A  general  algorithm  for  the  elimination  of  X  may  be 
written  as  n 

x»  •  'C'/C1'  -  T}WW*,  /•  *  * 1 . N  <»-»» 


A*..  a*-K  a'^Ma^Va""1) 
x/  x/  xn  n/  7  nit 


i,j  •  n  ♦  ( 4-21) 


8"  •  8”'1-  A”-1 (8rt"l/An‘l) 
x.  x  x.«  n  «H 


x  »  «  ♦  J , . ,N  (4-22) 


Equations  4-20,  4-21,  and  4-22  may  be  rewritten  in  compact 
form: 


■  *.  -  JVy  *  • "  *  ' . N 

x,  /  •  it  ♦  ? , . . .  . ,  N 


a1?...  a?:1  -  a1?:1  c 


M  x/ 


xn  it/ 


8.-8.  *  A .  1  0 

X.  X  XK  M 


x  •  n  ♦  i I W 


(4-23) 

(4-24a) 

(4-24b) 


,4-12 


where 


Pn  -  «;-/*- 1  C-25a) 

Cny  ■  <»-*»> 


After  the  above  procedure  is  applied  W- I  times  the  original 
set  of  equations  is  reduced  to  the  following  single  equa¬ 
tion 


8 


N-I 

W 


which  is  solved  directly  for  X^ 


XW 


jM-l 

N 


/A 


N-  / 
VN 


In  terms  of  the  previous  notation,  this  is 


XN  »  (4-26) 

The  remaining  unknowns  are  determined  in  reverse  order 
by  the  repeated  application  of  equation  4-23. 

4-5,2  Simplification  for  Band  Matrices 

The  stiffness  matrix  was  placed  in  a  Hband"  form  which 
resulted  in  the  concentration  of  the  elements  of  the  stiff¬ 
ness  matrix  along  the  main  diagonal.  Therefore,  the  fol¬ 
lowing  simplifications  in  the  general  algorithm  (equations 
4-23,  4-24a,  and  4-24b)  were  possible: 


x«  ■ 

P»  - 

\  -X. 

i  • 

it  ♦ 

♦  M-  / 

(4-27) 

a"  .  - 

a”'1  c  . 

-tit  n/ 

i.s  • 

n  ♦ 

♦  il-J 

(4-28) 

nit 

8  i  ' 

ntt-l 

Bi  - 

A1?’1  0 
*4 .«  n 

i  • 

it  ♦ 

^  «  «  *  »  i  H 

♦  W-  1 

(4-29) 

4-13 


♦ 


where  M  is  the  band  width  of  the  matrix. 

The  number  of  numerical  operations  can  further  be 
reduced  by  recognizing  that  the  reduced  matrix  at  any 
stage  of  procedure  is  symmetric.  Accordingly,  equation 
4-28  may  be  replaced  by  the  following  equation: 


i  •  n  *  f ,....,  a  ♦  M- I 
j  •  l, . . . . ,n  ♦  M- I 


(4-30) 


since 


The  number  of  numerical  operations  required  for  the 
solution  of  a  band  matrix  is  proportional  to  MM2  as  com¬ 
pared  to  M3  which  is  required  for  the  solution  of  a  full 
matrix.  Also,  the  computer  storage  required  by  the  band 
matrix  procedure  is  MM  as  compared  to  M2  required  by  a 
set  of  M  arbitrary  equations. 

4-5.3  Stresses 

Once  the  nodal  displacements  are  known,  it  is  pos¬ 
sible  to  determine  the  stresses  at  these  points.  For 
design  purposes  the  following  forces  are  required: 

Mx  ■  Moment  per  unit  length  in  the  x  direction 
-  Moment  per  unit  length  in  the  y  direction 
•  Twisting  moment  at  the  point 

Q  ■  Shear  force  per  unit  length  in  the  x  direc¬ 
tion 

Q  •  Shear  force  per  unit  length  in  the  y  direc- 
9  tion 

ax  *  Inplane  stress  in  the  x  direction 
■  Inplane  stress  in  the  y  direction 
t x^  ■  Inplane  shear  stress  at  the  point 
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The  moments  are  obtained  from  Equation  4-31 


where  P 


-p 

-up 

0 

U'XX  | 

-Up 

-P 

0 

< 

^•yy 

0 

0 

(l-u)P 

^,*y 

U} 

I  2 ( I -u  2 ) 


(4-31) 


The  shear  forces  are  obtained  from  Equation  4-32. 


>  -  -0 

w 

f 

J 

(4-32) 


Since  the  displacement  w  at  any  point  is  represented 
by  the  12-terra  polynomial  (Equation  2-l),  we  can  say 

«>(*,!/)  «  [  I  |x  J.'f  |  x2 1  x?/ Jtf2  |x3  j  x2'/|  xi/2  |x3  j  x3ff  1  xr/ 3  j 

or  {«;)  »  Cd][N](4)  (4-33) 

The  matrix  [N]**1  is  S3  shown  in  Figure  3-5, 


From  Equation  4-31,  4-32  and  4-33 


V 

1  u  0  0  0 

My 

u  10  0  0 

<?x 

1 

• 

0  0  o-f  0  0 

0  0  0  1  0 

r 

•c 

_ 

0  0  0  0  1 

,  we  get 

d'XK 

d,yy 

3/3x?2<* 


MU) 


(4-34) 
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Multiplying  the  derivatives  of  the  second  matrix 
by  the  first  and  assuming  u  •  0.2 ,  we  get 


fMx 

<  “u 

I8*' 


0 

0 

0 

1. 

0  .4 

5x  Ztj  ,4x 

1.2y 

6xy 

>.*»/] 

0 

0 

0 

.4 

0 

2 

1 . 2  x  4  y  2  x 

6y 

1. 2xy  6xy 

ff  “ 

0 

0 

0 

0 

1.6  0 

0  -3.2x-3.2y 

0 

-4.tx2-4.$y2 

0 

0 

0 

0 

0 

0 

6  0  2 

0 

by 

by 

0 

0 

0 

0 

0 

0 

0  2  0  . 

6 

bx 

bx 

cjux-JiL) 


For  the  normal  and  shearing  stresses,  we  have 


8 

a 

0 


ol 

l 

6 


vy 


®  X 
c 


[KHc} 


(4-35) 


where  a 

E 

uE  . 

n  * 

E  . 

lNo2)‘ 

s  * 

(  f-u 2)  ’ 

2 ( l*ul  * 

'out  ft,  ’ 

1 

\-b*y) 

ib-y) 

-y 

y 

0 

0  0 

0 

>  * 

0 

0 

0 

0 

l-a'x) 

-x  (a-x) 

X 

l"1*] 

i 

{-a+x} 

- X 

U-x)  x 

i-b*y) 

( b-y)  -y 

y 

or  < c )  -  [/]{u> 


(4-36) 
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Prom  Equations  4-35  and  4-36,  we  get 

{0}  *  [K][/]{u)  (4-37) 

These  moments,  shears  and  stresses  are  in  the  glo¬ 
bal  system.  Therefore,  for  proper  identification,  it  is 
necessary  to  convert  them  to  the  local  coordinates  of 
the  plates.  The  coordinate  transformation  matrix  is 
shown  in  Table  4-2. 


GLOBAL 

LOCAL  SYSTEM 

I 

Type  1 

Type  2 

Typ  3 

1 

1 

9 

5 

2 

-5 

-5 

*9 

3 

-9 

1 

1 

4 

11 

0 

10 

5 

10 

10 

0 

6 

0* 

11 

11 

7 

2 

12 

-6 

8 

-6 

-6 

-12 

9 

-12 

2 

2 

10 

14 

0 

13 

11 

13 

13 

0 

12 

0 

14 

14 

13 

3 

15 

-7 

14 

-7 

-7 

-15 

15 

-15 

3 

3 

16 

17 

0 

16 

17 

16 

16 

0 

18 

0 

17 

17 

19 

4 

18 

-8 

20 

-8 

-8 

-18 

21 

-18 

4 

4 

22 

20 

0 

19 

23 

19 

19 

0 

24 

0 

20 

20 

•This  and  all  zero3  appearing  here  are 
computer  controls  Ju3t  to  avoid  the 
local  coordinate  transformation. 


Table  4-2 

COORDINATE  TRANSFORMATION  COEFFICIENTS 
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4-6  Response  Spectra  For  Sheltered  Equipment 
4-6,1  General 

Equipment  sheltered  within  the  structure  will  be 
subjected  to  both  vertical  and  horizontal  movement  at 
its  point  of  attachment.  If  these  movements  cause  se¬ 
rious  damage  to  the  equipment,  it  may  matter  little 
that  the  structure  successfully  withstood  the  blast 
loading. 

Equipment  and  utilities  may  be  attached  to  the 
structure  by  means  of  a  flexible  support  system  (Fig-  ' 
ure  4-3), 


Steel 

Channels 

- — i  • 

Equipment 

I 

- 1 

a 

— ^  - 

L _ ! 

r 

m: 

-  , 

•  •*!  i 

J 

t 

i Igure  4-3  EXAMPLE  OF  A 
FLEXIBLE  HORIZONTAL 
MOUNTING 

Often,  this  is  done  in  order  to  isolate  the  item  from 
large  accelerations  or  may  result  simply  from  convention¬ 
al  attachne. procedures.  Because  the  mass-spring  system, 
coii.oosed  of  "‘quinnent  and  mounting,  can  vibrate  relative 
to  its  attachment  point,  its  peak  displacement  and  accel¬ 
eration  may  differ  appreciably  from  that  of  the  support¬ 
ing  structure.  por  design  purposes,  it  is  desirable  to 
have  a  graph  of  maximum  responses  at  an  attachment  point, 
for  a  wide  ran re  of  mount  stiffness-equipment  weight 
ratios.  This  graph  13  called  a  response  spectrum. 

A  socctrum  requires  that  the  reaction,  delivered  to 
the  supporting  structure  through  the  mounting,  does  not 
alter  the  timewise  movement  of  the  support  point.  Thus, 
the  weight  of  the  equipment  must  be  3rnll  compared  to  the 
supporting  member,  usually  less  than  10  percent. 
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4-6. 2  Definition  of  the  Response  Spectrum 

The  sinnle-derxee-of-freedom  system  considered  con¬ 
sists  of  a  rip;ld  mass,  m,  connected  to  the  point  of  at¬ 
tachment  by  both  a  weightless  spring  and  a  dashpot  (Fir,- 

ure  4-4). 


Ficure  4-4  SIMPLE 
ELASTIC  SYSTEM 


The  spring  and  dashpot  exert  restraining  forces  propor¬ 
tional,  respectively ,  to  the  relative  displacement  a.nc. 
relative  velocity  of  the  mass  and  attachment.  The  dis¬ 
placement  of  the  attachment  is  denoted  <j  and  the  spring 
deformation  by  x. 

The  equation  of  motion  for  the  system  (Figure  4-4) 
may  be  written, 

maj  +  c  i  +  kx  •  0  (4-33) 

Noting  that  w  •  y  *  i,  Equation  4-33  is  rewritten, 

mx  +  cx  +  kx  c  -mi/  (4-39) 

Or,  rewriting  in  terms  of  the  damping  ratio, 

mx  ♦  ZXumi*  kx  ■  -my  (4-40) 

For  a  specific  excitation,  mV,  of  a  simole  system  with 
a  particular  nerc'ntaqe  of  critical  dorpinr,  X,  the 
maximum  s-prinr  displacement,  x,  is  a  function  only  of  the 
circular  frequency  of  vibration  of  the  system, 
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A  plot  of  the  maximum  simple  system  response, 
whether  relative  displacement  x,  absolute  displacement 

w,  acceleration  :v,  or  spring  force  fcx,  against  the 
circular  frequency  of  vibration  w,  or  the  natural  fre¬ 
quency  of  vibration  /,  is  a  response  spectrum.  The 
most  useful  spectrum  is  that  of  relative  displacement , 

x.  In  the  following,  the  spectral  value  of  displace¬ 
ment  relative  to  the  point  of  attachment  will  be  desig¬ 
nated  as  (f. 

Maximum  values  of  spectral  pseudo-velocity,  l*,  and 
pseudo-acceleration,  A,  may  be  conveniently  defined  as, 

V  »  ui(l  =  UfU  (lJ-JJla) 

A  =  u2U  =  4*2f2U  (*U*llb) 

These  quantities  are  alternate  measures  of  the  maximum 
spring  deformation  U.  If  one  is  known,  the  other  two  are 
calculated  easily.  When  the  deformation  spectrum  is  plot¬ 
ted  on  a  logarithmic  scale  as  shown  in  Figure  *1-17,  all 
three  quantities  are  represented.  (See  Page  *1-36) 

Figure  *1-17  shows  the  response  spectra  for  one  of 
the  joints  in  the  sample  structure.  The  time-dependent 
acceleration  of  this  joint  is  calculated  by  the  pro¬ 
cedure  outlined  in  Section  *l-*l.  This  acceleration  be¬ 
comes  the  forcing  function,  !J,  of  Equation  *i-*iO.  The 
calculated  response  spectrum  is  simply  the  maximum 
spring  displacement  x  of  the  simple  clastic  system  of 
Figure  *!-*!,  The  general  characteristics  of  the  spectra 
arc  as  follows: 

1.  The  spectra  show  marked  peaks  at  a  fre¬ 
quency  of  20  cps,  corresponding  to  the 
predominant  box-structure  frequency  of 
response . 

2.  Peak  values  decrease  as  the  damping 
increases . 

3.  For  low  frequencies,  the  maximum  spectral 
displacements  (equipment  rattle-space) 
approach  the  peak  attachment  point  dis¬ 
placement  . 
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^4— 7  Sample  Dynamic  Analysis 


4-7.1  Orientation  of  the  Structure  and  Loading 

In  order  to  analyze  a  structure  in  the  most  ef¬ 
ficient  manner  with  respect  to  computer  time  it  is  neces¬ 
sary  to  orient  the  structure  in  such  a  way  that  the  num¬ 
ber  of  plates  in  the  x  direction  (see  Figure  2-1)  is 
equal  to  or  less  than  the  tj  plates  and,  the  number  of  <j 
plates  is  equal  to  or  less  than  tue  z  plates.  The  joint 
coordinates  are  Cartesian  and  follow  the  right-hand  rule. 

When  the  structure  has  been  oriented  to  satisfy  the 
mechanics  of  solution,  it  then  becomes  necessary  to  orient 
the  blast  pressures  and  ground  motion  to  be  compatible 
with  the  concept  of  the  problem.  There  are  six  possible 
cases  of  orientation  of  the  loading  which  are  numbered 
in  accordance  with  Figure  4-5.  The  Figure  depicts  the 
six  cases,  each  of  which  is  a  combination  of  joint  co¬ 
ordinates  and  load  coordinates.  All  permissible  condi¬ 
tions  can  be  satisfied  with  one  of  the  six  cases.  In  the 
case  of  a  model  not  oriented  In  its  natural  position  for 
mathematical  reasons,  there  may  be  more  than  one  way  in 
which  it  can  be  turned  but  at  least  one  way  will  be  sat¬ 
isfied  by  one  of  the  cases  of  Figure  4-5. 

4-7.2  Numbering  Sequence 

Each  plate  of  the  model  is  divided  into  four  rec¬ 
tangular  elements.  The  division  should  be  visualized  as 
being  accomplished  by  cutting  planes  entirely  through  the 
model.  Wherever  ^hree  planes  (plates  or  cutting  planes) 
intersect  there  is  a  joint.  Where  all  three  planes  are 
cutting  planes,  the  joint  is  a  dummy  joint.  Each  Joint 
is  given  a  number  by  the  computer.  These  numbers  are  in 
sequence  without  regard  to  whether  the  joint  is  actual 
or  dummy. 

The  joints  are  numbered  in  x,  y ,  z  sequence.  (Fig¬ 
ure  2-1).  The  computer  starts  on  the  Xf /  plane  with  the 
smallest  value  of  z.  The  numbering  is  row-wise  sequenced 
downward  and  continuing  through  the  z  planes  in  ascend¬ 
ing  order. 
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The  elements  are  numbered  in  a  manner  similar  to  the 
Joints.  In  each  plane  they  are  numbered  left  to  right, 
top  to  bottom.  The  numbering  starts  with  the  elements  in 
the  xj  plane,  starting  at  the  smallest  value  of  z,  and 
working  to  the  largest  value  of  z.  It  then  continues  in 
the  tjz  plane  starting  with  the  smallest  value  of  x  and 
working  through  to  the  largest  value  of  x.  Finally,  it 
takes  the  zx  planes  from  the  lowest  to  the  highest  value 
of  ij.  No  differentiation  is  made  between  the  real  and 
dummy  elements.  The  exploded  box  (Figure  4-6)  demon¬ 
strates  the  numbering  sequence. 

The  order  of  recognition  of  the  plates,  e.g.,  for 
inputting  the  plate  thicknesses,  is  the  same  as  for  the 
elements,  however,  the  dummy  plates  which  are  introduced 
in  the  analysis  are  ignored. 

4-7.3  Listing  of  the  Data  on  the  Input  Forms 

This  Sub-section  should  be  read  in  conjunction  with 
input  sheets  1  and  2  (Figures  4-9  and  4-io). 

Job  identification  is  an  alphameric  heading  occupy¬ 
ing  not  more  than  72  spaces. 

In  the  second  line  nx,ny  and  itz  are  fixed  point 
digital  entries  equal  to  the  number  of  cross  walls  in  the 
x,ij  and  z  directions  respectively.  If  a  calculation  of 
response  spectra  Is  required,  a  digit  is  placed  under  the 
heading  "Response  Spectra."  If  this  is  left  blank,  no 
spectra  will  be  calculated  and  input  sheet  No.  2  will  not 
be  required.  The  Case  Number  refers  to  the  orientation 
of  the  structure  and  loads  as  described  in  sub-section 
4-7.1  and  Figure  4-5. 

The  following  three  lines  are  for  the  lengths  of  the 
plates  in  the  x,tj  and  z  directions.  In  the  following 
lines,  the  thicknesses  of  the  plates  are  entered  in  se¬ 
quence  as  described  in  Sub-Section  4-7.2. 

In  line  9  are  entered  the  modulus  of  elasticity  and 
density  of  the  material  of  the  structure,  and  the  density 
and  thickness  of  the  soil  layer. 

The  next  two  lines  contain  the  data  required  to  de¬ 
fine  the  loads  due  to  the  blast  and  ground  motion,  which 
are  explained  in  Appendix  A. 
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In  line  12  are  entered  the  viscous  damping  factors 
and  the  numerical  integration  controls,  i.e.,  the  time 
interval  At  and  the  number  of  time  intervals. 

The  next  four  lines  are  used  to  identify  the  joints, 
if  any,  at  which  time  dependent  print  out  of  the  displace¬ 
ments,  velocities  and  accelerations  are  required,  and 
also  which  coordinate  direction  these  responses  are  re¬ 
quired.  The  coordinates  are  indicated  numerically,  i.e., 
x  »  1 ,  ij  *  2 ,  z  *  3.  Because  of  the  large  amount  of 
printed  data  that  this  computer  program  is  capable  of 
producing,  it  was  decided  to  limit  the  time  dependent  p-rint 
to  three  nodal  points  per  run. 

If  response  spectra  are  required,  input  sheet  No.  2 
must  be  used.  Spectral  velocities  can  be  obtained  for  any 
twenty  nodal  point  translational  displacements.  Four 
different  damning  factors  may  be  used.  The  frequencies 
at  which  the  spectral  velocities  are  required  are  listed 
row-wise  in  sequence  to  a  maximum  of  twenty-five. 

4-7. 4  Sample  Analysis 

The  modeling  of  a  structure,  numbering  of  elements, 
reorientation  of  loading  coordinates  and  entry  of  data  on¬ 
to  the  input  forms  is  demonstrated  in  the  sample  analysis. 
The  structure  is  shown  in  its  natural  position  in  Fig¬ 
ure  4-7  And  reoriented  with  reoriented  loading  in  Fig¬ 
ure  4-8,  The  numbering  of  the  elements  is  identical  with 
the  sample  of  Figure  4-6.  The  data  input  forms  are  shown 
In  Figures  4-9  and  4-10.  A  printing  of  the  input  is  shown 
In  Figure  4-11,  A  sample  page  from  the  output  of  load 
deflection  history  for  two  points  is  shown  In  Figure  4-12. 
Figures  4-13,  4-14  and  4-15  show  the  maximum  displacements, 
velocities  and  acclerations .  After  each  of  these  list¬ 
ings,  the  time  at  which  these  maxima  occured  is  also  print¬ 
ed,  Fi  ure  4-l6  shows  the  maximum  positive  stresses  in 
the  elements  in  local  coordinates.  The  complete  printout 
lists  the  moments  M  ,  M  and  M  ,  the  forces  Q  and  Q 
and  the  stresses  o  a  and  TX(n.  This  listing  is  follow¬ 
ed  by  a  listing  of*the^tlmes  at  hvhich  these  maxima  oc¬ 
cured,  a  listing  of  the  maximum  negative  stresses  and  the 
times  at  which  the  maximum  negative  stresses  occured. 

Figure  4-17  shows  response  spectra  for  the  midpoint  of 
the  base  slab  in  the  vertical  direction. 
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MAXIMUM  ACCELERATIONS 


MAXIMUM  POiiTlVt:  STRESSES  I n  LOC  AL  COOPOINA TES 


ELEMENT  *  TV ?t  »  JOINT  #  "  "MOM>'MT-X . MO'iPNT-Y 


a  I 

b  i 

c  l 

d  1 

e  2 

f  2 

fS  ‘  2 

h  2 

I  3 

J  3 

k  3 

1  3 

m  4 

n  4 

o  4 

P  4 


QX 


a  O.OOOO+OOO 
b  0 .0000*000 
C  5.9934+000 
d  5,9934+000 
e  1.9474+001 
f  1.9474+001 
K  1,1378-002 
h  1. 1378-002 
1  7,0713+000 
J  7,6713+000 
k  7, ?9o 1 ♦ QOO 
1  7 , ?9o 1  +  OOO 
m  2,9172-002 
n  2,9172-002 
o  5,9292+000 
P  5,9292+000 


1  4 

1  S 

1  1 

l  2 

1  5 

1  6 

1  2 

1  3 

1  7 

1  8 


1  9 

l  5 

1  b 
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0.0000+000 
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2 . 01 1  5  ♦  n  0  1 
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l , 652^-00  1 
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0,0000+000 

0,0000+000 

0.00O0+040 

6.7959-002 
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MAXIMUM  POSITIVE  STRESSES 


Frequency  (cps) 

Figure  4-17 

RESPONSE  SPECTRA  FOR  JOINT  41  (z  Direction) 
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Section  5 


DYNAMIC  RESPONSE  ANALYSIS  -  MOM-LINEAR 
5-1  General  Procedure 

The  non-linear  dynamic  analysis  is  carried  out  by 
a  step-by-step  procedure  in  which  the  structure  is  as¬ 
sumed  to  respond  linearly  during  each  short  time  in¬ 
terval.  The  elastic  properties  may  be  changed,  however 
from  one  Interval  to  the  next;  thus,  the  non-linear 
response  is  obtained  as  a  sequence  of  linear  responses 
of  successively  differing  systems. 

The  analysis  procedure  involves  the  repeated  appli 
cation  of  the  following  steps  for  each  successive  time 
interval: 

First:  the  stiffness  matrix  appropriate  to  the  time 

interval  is  evaluated,  based  on  the  forces 
existing  in  the  plate  elements  at  the  begin¬ 
ning  of  the  time  interval. 

Second:  changes  in  displacements  of  the  elastic  struc 

ture  are  computed,  assuming  the  accelerations 
to  vary  linearly  during  the  interval. 

Third:  these  incremental  displacements  are  added  to 

the  deformation  state  existing  at  the  begin¬ 
ning  of  the  interval,  to  obtain  total  element 
deformations . 

Finally:  based  on  these  element  deformations,  plate 

forces  are  computed  from  which  the  stiffness 
matrix  appropriate  to  the  next  time  interval 
may  be  determined. 

5-2  Failure  Criterion 
5-2.1  Method  of  Failure 

There  is  no  well  established  theory  of  failure  of 
plates  under  the  combined  effect  of  bending  and  inplane 
forces.  Therefore,  in  this  study  an  assumption  was  made 
that  the  plates  would  fail  in  bending  first,  and  once 
having  failed,  a  plate  could  not  recover  its  bending 
stiffness.  Also,  a  plate  is  assumed  to  continue  to 
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contribute  elastically  to  the  overall  stiffness  of  the 
structure  through  its  inplane  stiffness  after  losing 
its  bending  stiffness. 

Each  element  can  lose  its  bending  stiffness  in¬ 
dependent  of  any  other  element.  The  whole  structure  is 
considered  to  have  failed  when  the  application  of  load 
is  such  that  a  sufficient  numher  of  elements  have  failed, 
so  that  the  stiffness  of  the  structure  becomes  singular. 

The  stiffness  matrix  will  become  singular  when  all 
the  elements  coinciding  at  a  node  point  will  not  provide 
additional  flexural  constraint  thereby  permitting  infi¬ 
nite  displacements  at  that  node. 

5-2.2  Failure  Assumption 

Figure  5-1  shows  the  Johansen's  yield  line  theory11 
applied  to  a  uniformly  loaded  rectangular  plate  clamped 
on  all  four  edges. 


figure  5-1 
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Based  unon  vield  line  theory,  a  kinematically  com¬ 
patible  system  of  forces  gives  an  upper  bound  to  the 
solution  whereas  using  any  statically  admissible  system 
will  give  a  lower  bound.  Investigations  have  shown 
that  the  difference  between  the  two  methods  of  solution 
is  small,  and  that  the  yield  line  theory  gives  conserv¬ 
ative  answers.  (See  Section  6) 

The  angle  $  is  determined  so  that  a  minimum  value 
of  load  p  is  obtained  which  will  cause  the  plate  to 
reach  its  yield  canacity  M  , 

*  p 

Assuming  the  moment  capacity  is  the  same  for  pos¬ 
itive  ann  negative  moments  in  both  the  X  ana  V  directions 
it  can  be  said, 

(1)  External  Work 

fa 2  fa 2  fa  I 

t»'e  =  pj  tan$  +  pj  tan$  +  2p  j  (a-fatan^)  j 

=  4’tan*  ♦  4  U-Han*>  (5_1> 


b-3 


JS'I  i!  ,  <*  /  I 

p  3tan<j>  7btan2$  Ttand 


r .  JI 

S’l  .  i  .  +  r  =  0 

p  t  m  ij  b: 


or  tan<j> 


-  s.  ♦  ♦  j 


(5-4) 


5-2.3  Failure  Condition 

Since  the  finite  elements  are  equilibrated  to  nodal 
points  at  the  four  corners,  it  is  necessary  to  make 
some  assumption  as  to  the  development  oh  the  yield 
moment.  Che  method  elected  was  to  average  the  moments 
ir.  the  $  direction  at  all  four  corners  and  compare  the 
average  to  moment  capacity  of'  the  element. 

Using  ilohr's  circle,  the  moment  at  a  point  vrhich 
makes  an  angle  <j>  with  the  V  axis  is , 

=  Mycos2<}>  -  ( x.js in2 <J>  +  ’lxsin2$  (5-5) 

and  the  failure  condition  is, 

■V  <  7  +  V  +  \k  *  \l)  (5'6) 


5-3  Sanole  Non-T, inear  Analysis 

A  single  cell  box  structure  was  defined  as  shown  in 
^ l mure  5-2, 

The  structure  was  40  hf»ct  square  by  20  reet  deep 
with  walls  0,3  reet  thick,  Different  moment  cnnacitlea 
were  assigned  to  each  elate  so  that  the  e^ect  of  suc¬ 
cessive  "allures  could  bo  demonstrated. 

Figure  5-3  shows  the  input  form  and  Figure  5-4 
shows  a  printing  of  the  input. 
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rM.«*uro  5-2 

sa'tlf  nr  t^otid  .t^ucturf: 


27 

\ 

> 


Figures  5-5,  5-6,  and  5-7  are  selected  sheets  from 
the  output,  Figure  5-5  shows  the  Load  Deflection  History, 
After  each  lncrcnent  of  tine  the  moments  In  each  plate 
are  compared  to  the  assumed  average  moment  capacity  of 
each  plate.  When  any  plate  has  failed,  the  plate  number 
is  printed  on  the  output.  The  process  continues  until  a 
sufficient  number  of  plates  have  failed  to  produce  an 
unstable  structure,  at  which  point  the  computer  halts  any 
further  analysis  and  prints  out  maximum  displacements 
(Figure  5-6),  velocities  and  accelerations.  The  maximum 
positive  moments  and  maximum  normative  moments  (Figure  5-7) 
arc  also  printed  out.  The  times  at  which  those  maxima 
occur  are  also  printed. 
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”l~urc  5-3 
lO'i-LCN'I'AR  nWRAM  - 


.'it .  CL5/c/'_). . 
150.  2, 


oynamic  non-linear  analysis  of  box  structures 
sample  PROBLEM  ...  SINGLE  box  DYNAMIC  NON-LINEAR  RESPONSE  ...  10/12/66 
NO#  OF  JOINTS 

X  DIRECTION  3  2 

Y  DIRECTION  a  2 

Z  DIRECTION  3  2 

MODULUS  OF  ELASTICITY 
STRUCTURE  »  432000.00  KSF 

SOIL  =1255444.72  KSF 
RO  =  150.00  LB/CU.FT. 

KO  ■  0.33 

LENGTHS  OF  PLATES 

X  DIRECTION  3  40,00 

Y  DIRECTION  a  20.00 

Z  DIRECTION  a  40.00 

plate  thicknesses 

0,30  0.30  0,30  0,30  0,30  0.30 

MOMENT  CAPACITY  <K»FT/FT) 

20,00  20.00  15.66  15.00  100.00  500.00 

blast  and  ground  motion  characteristics 


W(MT)  8 

20.00 

PSO(PSl)= 

15*00 

al^ha  a 

1.36 

BETA  a 

8.40 

A  ■ 

0*67 

DP* ISEC)= 

0*08 

U < KPS)  3 

15000.00 

R (FT)  3 

12500.00 

OELH (FT) 3 

B300.00 

VL(FPS) 3 

17000*00 

i/L^FPS)3 

11700.00 

TM/T2  a 

0,83 

SP.GRAV.3 

2*69 

P  FACTOR3 

4*00 

PHASING  3 

-0*000000 

damping  factors 

alpha  a  loo.oooooo  beta  * 


40  TIME  INCpEMENTS 

TIME  INCREMENT  a  0.0005 

6.000025 


Figure  5-^ 
PRINTING  OF  INPUT 
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LOACUOEFLECTION  HISTORY 


IOAO  IN  KIPS*  DEFLECTION  IN  FT*  TIME  IN  SEC 

5-2  ll-Y  1 3-X 


time 

LOAD  DEFLECTION 

load 

deflection 

load 

deflection 

0,00050 

5,3  0.000000 

155,9 

0.000011 

71.1 

0.000001 

N£W  FAILED 

PLATES0*  NONE 

• 

0,00100 

21,4  0.000000 

311,7 

0.000090 

140.6 

0.000005 

NEW  FAILED 

PLATES'40  NONE 

0.000014 

0.00150 

47,4  0.000002 

463.6 

0,000301 

136,6 

NEW  FAILED 

PLATES**  NONE 

0.00200 

72,6  0.000005 

607.9 

0,000703 

132,8 

0,000022 

NEW  failed 

PLATES00  NONE 

0.00250 

97,1  0.000009 

752,7 

0.001352 

129.1 

0,000024 

NEW  FAILED 

plates°°  none 

0.00300 

120.1  0.000013 

786.1 

0.002290 

125.6 

0.000019 

NEW  FAILED 

PLATES00  none 

0*00350 

130.3  .0.000016 

764,6 

0.003522 

122.2 

0,000011 

NEW  FAILED 

PLATES°°  NONE 

0.00400 

128.1  0.000018 

743.8 

0,005031 

119.0 

0.000006 

NEW  FAILED 

plates°°  none 

0.00450 

124,6  0*000018 

723.9 

0.006793 

115,8 

0.000007 

NEW  FAILED 

PLATES00 

15 

16 

0,00500 

121.2  0.000017 

704,7 

0.008790 

112.8 

0.000015 

NEW  FARED 

plates™ 

3 

11 

0.U0550 

U8.0  O.OOOOlfl 

686.3 

0.011003 

lo9.9 

0.000019 

NF.W  failed 

PLaTeS°° 

4 

12 

23 

24 

0,00600 

114,9  0.000017 

668,6 

0.013415 

107,2 

0.000017 

NEw  FAILEO 

PLATES™  NONE 

0.00650 

112.0  0.000012 

651 .6 

0.016012 

104,5 

0,000009 

nf.w  failed 

PLATES*0  none 

0,00700 

109,1  0.000005 

635,1 

0.018779 

101.9 

-0.000000 

new  failed 

PLATES00 

13 

14 

21 

22 

FaILEO  PLATES  PRODUCE  SINGULAR  STIFFNESS  MATRIX. 


STRUCTURE  BECOMES  UNSTABLE 


Figure  5-5 

LOAD-DEFLECTION  HISTORY 
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MAXIMUM  DISPLACEMENTS 
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MAXIMUM  DISPLACEMENTS 


MAXIMUM  NEGATIVE  STRESSES  IN  LOCAL  COORDINATES 


ELEMENT  •  TYPE  #  JOINT 

*  MOMENT-X 

MOMENT-Y 

moment-xy 

1 

1  4 

-5.5221*000 

-2,9776*001 

-2,7428-002 

l 

1  5 

-8,5711-001 

-2,2518*000 

-8,4298-003 

l 

1  1 

0. 0000*000 

0.0000*000 

-3.2826-001 

1 

1  2 

0,0000*000 

0,0000*000 

-3,4805-001 

2 

1  5 

—  1 , 0471-00 1 

-2,1077*000 

-3.7118*000 

2 

1  6 

-4,1761*000 

-1,9107*001 

-3,5308*000 

2 

1  2 

0,0000*000 

0,0000*000 

-5,8600-001 

2 

1  3 

-3,3463-003 

-1,7338-004 

-4. 0493-661 

3 

1  7 

-3,9758*000 

-1,9922*001 

-4,8388*000 

3 

1  8 

-8,0417*000 

-4,0164*001 

-4,7730*000 

3 

1  4 

0,0000*000 

0,0000*000 

0.0000*000 

3 

1  s 

-4,6388-002 

0,0000*000 

0.0000*000 

4 

1  8 

-1 ,2862*001 

-6,4319*001 

0,0000*000 

4 

1  9 

-2,2483*000 

-1,1246*001 

-4,7270-667 

4 

1  S 

0,0000*000 

0,0000*000 

-8,3182-001 

4 

1  6 

-6,2172-002 

-9,5414-004 

-9,8023-001 

5 

1  13 

0,0000*000 

0,0000*000 

0,0000*000 

5 

1  14 

0,0000*000 

0,0000*000 

0.0000*000 

5 

1  10 

0,0000*000 

0,0000*000 

0,0000*000 

5 

1  11 

0,0000*000 

0,0000*000 

0.0000*000 

6 

1  14 

0,0000*000 

0,0000*000 

0.0000*000 

6 

1  15 

0,0000*000 

0,0000*000 

o.oooo*6oo 

6 

1  11 

0,0000*000 

0,0000*000 

0,0000*000 

6 

1  12 

0,0000*000 

0,0000*000 

0,0000*000 

7 

1  16 

0,0000*000 

0,0000*000 

0.0000*000 

7 

1  17 

0,0000*000 

0,0000*000 

0.0000*000 

7 

1  13 

0,0000*000 

0,0000*000 

0,0000*000 

7 

1  14 

0.0000*000 

0,0000*000 

0.0000*000 

8 

1  17 

0,0000*000 

0.0000*000 

0,0000*000 

8 

1  18 

0,0000*000 

0,0000*000 

0,0000*000 

8 

1  14 

0,0000*000 

0,0000*000 

0,0000*000 

8 

1  15 

0.0000*000 

0,0000*000 

0.0000*000 

9 

1  22 

-2,5286-002 

0,0000*000 

-4.6697*000 

9 

1  23 

0,0000*000 

0,0000*000 

-4,5583*000 

9 

l  19 

-5,6850*000 

.2,8536*001 

-1,3760*000 

9 

1  20 

-1,6973*000 

-8,3749*000 

-1,2646*000 

10 

1  23 

-1,8730-001 

0,0000*000 

-7,9570-004 

10 

l  24 

0.0000*000 

0.0000*000 

-2,3633-002 

10 

1  20 

-1.6542*000 

-8.3662*000 

-5,3776*001 

10 

1  21 

-3.6424*000 

-1*8131*001 

-6*7632-001 

Figure  5-7 

MAXIMUM  NEGATIVE  STRESSES 
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Section  6 


SYNOPTIC  REVIEW  OF  SELECTED  LITERATURE  PERTAINING  TO  THE 
FAILURE  MODES  AND  FAILURE  LOADINGS  OF  TWO-WAY  SLABS 

6-1  Introduction 

■  In  the  search  for,  and  review  of  publications,  emphasis 
was  placed  on  reports  of  the  results  of  load  tests  to 
failure  and  comparisons  of  test  results  with  predict¬ 
ions  calculated  by  analytical  methods.  The  region  of 
interest,  in  the  load-deflection  history  of  slabs, 
starts  at  the  beginning  of  yield  in  the  reinforcement 
and  ends  at  collapse  of  the  slabs.  Slabs  are  assumed 
to  have  less  than  balanced  reinforcement. 

All  of  the  tests,  of  which  reports  were  examined,  were 
made  with  static  loads.  Where  the  location  of  the  test 
was  reported,  it  is  given  in  this  review.  Generally, 
locations  are  not  reported  but  inasmuch  as  the  tests 
were  made  in  laboratories,  it  is  assumed  that  they  were 
facilities  of  the  activity  with  which  the  author  was 
connected. 

In  the  listing  of  references,  those  publications  which 
were  considered  to  be  significant  contributors  of  know¬ 
ledge  in  the  region  of  interest  are  abstracted.  Many 
other  reports  were  examined  and  those  addressed  to  ob¬ 
jectives  which  bear  to  some  lesser  extent  on  the  region 
of  interest  are  listed  but  not  abstracted. 

6-2  Synoptic  Review  of  Selected  Literature 

The  pertinent  phenomena  reported  In  the  region  of 
interest  are: 

1)  The  yielding  of  the  reinforcing  steel. 

2)  The  arching  action  of  the  slab. 

3)  The  improvement  of  the  plastic  deformability 
of  the  concrete  in  the  triaxially  compressed 
zone. 

4)  The  mobilization  of  tensile  membrane  forces 
In  the  suspended  net  of  steel. 

The  excitation  of  interest  in  the  post-yield  history 
of  overloaded  slabs  appears  to  have  been  triggered  by 
Ocklestons 1 "report  of  his  load  tests  to  destruction  of 

^References  designated  by  superscript  in  this  Section  6 
are  found  on  pages  6-6  to  6-12. 


some  two-way  slab  floor  panels  in  a  reinforced  concrete 
building  in  Johannesburg,  South  Africa,  which  was  avail¬ 
able  because  it  was  scheduled  for  demolition. ' Ockleston 
reported  cracking  patterns  agreeing  with  the  Johansen 
Yield  Line  Theory,  but  that  he  obtained  failure  loads 
over  2h  times  analytical  predictions  according  to  the 
theory.  Approximately  2\  years  later  (June  1958)  he 
reported2  his  re-examination  of  his  test  results  and 
his  success,  with  a  simple  theory  which  he  produced, 
in  rationalizing  the  high  failure  loadings  as  result¬ 
ing  from  the  superposition  of  plastic  hinge  moments  and 
the  resisting  moment  produced  by  arching  action  of 
compressive  axial  force  induced  in  the  slab  by  large 
deflections  associated  with  the  functioning  of  plastic 
hinges . 

A  somewhat  idealized  version  of  arching  action  in  a 
slab  is  depicted  in  Figure  6-1.  It  is  assumed  that 
the  slab  is  under-reinforced,  as  is  normally  true, 
so  that  the  concrete  behaves  elastically  above  the 
yield  moment  of  the  steel. 

Relatively  large  deflection  of  the  slab  after  the  form¬ 
ation  of  plastic  hinges  imposes  compressive  strains  for 
the  simple  reason  that  the  diagonal  distance,  d ,  between 
compressive  stress  block  centroids,  is  greater  than 
the  horizontal  distance,  1/2. 

It  is  readily  apparent • that  these  compressive  strains 
will  not  be  imposed  unless  the  structure  surrounding 
the  yielding  slab  is  stiff  enough  to  function  as  an 
effective  abutment.  It  is  also  apparent  that  the 
limit  of  arching  action  is  reached  when  the  lever  arm 
reaches  zero  and  that,  with  small  additional  deflection, 
the  lever  arm  becomes  negative  and  the  stored  inplane 
strain  energy  performs  work  contributing  to  plastic 
hinge  rotation.  This  produces  additional  deflection 
with  rapidly  diminishing  load,  transforming  work  already 
done  into  tensile  membrane  strain  energy  in  the  steel, 
which  is  already  past  its  yield  point.  With  sustained 
load,  the  rupture  stress  in  the  steel  is  likely  to  be 
reached,  particularly  with  large  depth  span  ratios  and 
very  stiff  abutting  structure.  With  impulse  loads, 
the  work  of  which  is  quickly  transformed  to  kinetic 
energy,  the  kinetic  energy  remaining  when  the  steel 
becomes  the  sole  survivor  may  be  within  the  capacity  of 
the  steel  to  accept  it  below  the  rupture  stress.  The 
worth  of  the  suspended  steel  net  depends  considerably 
on  the  stress-strain  properties  of  the  steel  above  its 
yield  point  stress. 
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Figure  6-1 

ARCHING  ACTION  IN  A  YIELDING  SLAB 


The  characteristic  load  deflection  curve  is  shown  in 
Figure  6-2. 


Figure  6-2 


CHARACTERISTIC  LOAD  DEFLECTION  CURVE  OF  REINFORCED  SLAB 
WITH  ARCHING  ACTION  MOBILIZED  BY  LATERAL  BOUNDARY 

STIFFNESS 
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The  foregoing  synopsis  of  load-deflection  history  is  a 
syllabus,  to  a  large  extent,  of  the  findings  from  ex¬ 
amination  of  References  1-12,  inclusive.  Tests  re¬ 
ported  by  Park3  and  by  Christiansen4  manifested  the 
existence  of  arching  action  and  emphasized  its  de¬ 
pendence  on  lateral  boundary  stiffness.  Christiansen* s 
results  with  two  quite  widely  differing  cube  strengths 
pointed  up  the  sensitivity  of  arching  action  to  the 
inplane  stress-strain  relation.  Self5  reported  tests 
with  cold-rolled  and  hot-rolled  steel  which  suggested 
the  effect  of  strain  hardening  on  the  plastic  hinge 
resisting  moment.  Gamble  and  Sozen7  reviewed  the  lit¬ 
erature  and  studied  the  results  of  five  scale  model 
tests  made  at  the  University  of  Illinois  by  Seiss7,  and 
advanced  five  reasons  for  the  calculation  of  failure 
loads  significantly  less  than  test  results.  These  are: 

1)  Incorrect  determination  of  steel  strength 

2)  Strain  hardening  of  steel 

3)  Erroneous  location  of  yield  line 

4)  Arching  Action 

5)  Effects  of  deformation  of  the  structure 

Undoubtedly,  the  increase  in  load  capacity  above  the 
yield  point  in  the  steel  is  a  superposition  of  strain 
hardening  and  arching  effects.  Hodge  and  Perrone6, 
Hillerborg8  9,  Crawford10  and  Kemp11  report  finding  by 
means  of  theory  that  the  yield  line  method  gives  an 
upper  bound  and  that  there  is  a  lower  bound  that  can 
be  determined  by  finding  a  distribution  of  moments 
which  satisfies  the  boundary  conditions  and  the  equi¬ 
librium  equation9.  Kemp  found  the  lower  bound  to  be 
only  slightly  less  than  the  upper  bound,  but  found  the 
yield  lines  to  be  an  infinite  set  of  curves.  These 
theoretical  findings  regarding  the  failure  load  are 
not  borne  out  by  the  tests  reported  upon,  in  which  test 
results  were  higher  than  calculated  values.  Gamble  and 
Sozen7  found  thi3  to  be  true  in  the  literature  they  re¬ 
viewed  and  the  results  of  scale  model  tests  they  stud¬ 
ied.  Relative  to  their  advancement  of  incorrect  de¬ 
termination  of  steel  strength  as  one  of  the  reasons  for 
low  calculated  failure  loads,  the  authors  of  the  publi¬ 
cations  which  have  been  examined,  in  which  analytically 
predicted  and  experimental  results  were  compared,  did 
not  report  whether  or  not  they  had  determined,  by  test, 
the  yield  point  of  the  specific  reinforcement  earmarked 
for  use  in  experimental  loadings.  The  yield  point 
stresses  specified,  according  to  steel  grades,  are  the 
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producers'  guaranteed  minima,  which  are  exceeded  approxi¬ 
mately  15  percent  by  the  mean  value. 

Brotchic,  Jacobson  and  Okubo12  report  the  results  of 
the  comprehensive  tests  and  investigation  made  at  the 
Massachusetts  Institute  of  Technology  under  contract 
with  the  U.  S.  Naval  Civil  Engineering  Laboratory, 

Port  Hueneme,  California.  This  recent  (August  1965) 
supported  work  in  the  field  appears  to  have  been  a 
well  integrated,  comprehensive  experimental  investi¬ 
gation  of  the  behavior  of  two  way  square  slabs,  sup¬ 
ported  around  the  four  edges,  with  and  without  moment 
restraint.  They  were  restrained  against  lateral  dis¬ 
placement,  but  with  strain  gauge  cells  which  permitted 
the  restraining  forces  to  be  measured.  The  report 
presents  straightforward  analytical  equations  which, 
with  a  few  exceptions,  agree  reasonably  well  with  the 
results  of  the  experiments.  A  feature  of  the  report, 
not  found  in  other  literature  examined,  deals  with  the 
effect  on  concrete  strength  of  the  triaxial  compression 
in  the  positive  moment  hinge  area,  comprised  of  axial 
forces  in  two  directions  and  the  vertical  force  of 
applied  load.  The  authors  conclude  that,  "The  com¬ 
bination  of  increased  load  capacity  and  Improved 
behavior  produced  by  arching  action  is  sufficiently  large 
and  sufficiently  predictable,  to  warrant  its  serious 
consideration  in  design."  They  suggest  Its  application 
to  the  roof  and  walls  of  reinforced  protective  structures 
of  various  kinds,  eg.  shelters,  block  houses,  etc., 
particularly  subterranean.  They  caution  that,  "tensile 
membrane  action  may  be  considered  only  in  situations  where 
excessive  deflections  are  tolerable  and  reuse  is  not 
necessary:  such  as,  protective  slabs  under  debris  loading." 

Regarding  the  predictability  of  the  increased  load  capacity 
produced  by  arching  action,  the  authors  did  not  address 
themselves  to  the  problem  of  determining  the  lateral 
stiffness  of  the  boundary  of  the  loaded  slab  in  an  exist¬ 
ing  structure.  Arching  action  is  sensitive  to  lateral 
stiffness  at  the  boundary  and  this  parameter  needs  to  be 
estimated  rather  closely.  A  structure  may  be  designed 
with  boundary  lateral  stiffness  provided  in  such  way  that 
its  effectiveness  may  be  reliably  verified  by  analytical 
methods,  but  such  design  may  not  be  found  in  existing 
structures. 

The  capability  of  the  tensile  membrane  to  provide  pro¬ 
tection  against  catastrophic  failure  under  destructive 
dynamic  loading  of  buried  structures  is  subject  to  certain 
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reservations  which  demand  thorough  consideration.  Tensile 
membranes  which  are  the  surviving  intact  components  of 
external  horizontal  slabs  will  be  subjected  to  the  sus¬ 
tained  gravity  load  of  earth  overburden  and  residual 
kinetic  forces.  Cracked  external  walls  with  relatively 
large  deflections  will  be  axially  loaded  with  the 
gravity  reaction  of  earth  overburden.  If  the  concrete 
is  not  too  badly  crushed,  this  force  will  provide  some 
resistance  against  the  lateral  forces.  If  the  concrete 
toward  the  outside  face  is  ineffective  in  compression,  the 
line  of  force  may  be  so  far  toward  the  inside  face  that 
the  gravity  reaction  will  contribute  to  collapse  of  the 
walls.  The  tensile  membrane  has  no  capability  to  resist 
collapse  from  such  cause.  The  lateral  force  will  be 
comprised  of  residual  kinetic  force  and  lateral  pressure 
of  disturbed  soil.  The  tensile  membrane  has  capacity 
within  its  strength  limitation  to  resist  these  forces, 
provided  the  vertical  force  does  not  cause  collapse. 

6-3  Selected  References  with  Abstracts 
of  More  Pertinent  Publications 

1.  Ockleston,  A.  J.,  Load  Tests  on  a  Three  Story 
Reinforced  Concrete  Building  in  Johannesburg ,  The 
Structural  Engineer  (London)  October  1955,  P .  30V.’ 

Full  scale  static  loading  tests  to  failure  were  run 
on  two  l6ft,by  13ft. 6in,  reinforced  concrete  slab  panels 
which  were  part  of  a  beam  and  slab  floor  system  27ft. 
wide  by  over  100ft.  long.  Slabs  were  bounded  by 
transverse  main  beams  spaced  l6ft.  apart  and  sec¬ 
ondary  longitudinal  beams  spaced  13ft.  6in.  apart. 
Cracking  pattern  agreed  with  Johansen  Yield  Line 
Theory  but  failure  loads  were  over  2h  times  the 
analytical  predictions. 

2.  Ockleston,  A.  J.  ,  Arching  Action  in  Reinforced  Con¬ 
crete  Slabs ,  The  Structural  Engineer  (London) 

June  1958,  P.197. 

Ockleston  re-examined  the  results  of  his  earlier 
tests  (ref.  1)  and  succeeded,  with  a  simple  theory 
he  produced,  in  rationalizing  the  high  failure 
loads  as  resulting  from  the  superposition  of  plastic 
hinge  resistance  and  the  resisting  moments  afforded 
by  compressive  membrane  forces. 
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3.  Park,  R.,  The  Lateral  Stiffness  and  Strength  Re¬ 
quired  to  Ensure  Membrane  Action  at  the  Ultimate 
toad  of  a  Reinforced  Concrete  Slab  and  Seam  Floor, 
Magazine  of  Concrete  Research  (London)  Vol.  17, 

No.  50,  March  1965,  PP. 29-35. 

Sand-cement  mortar  models  of  a  nine  panel  slab 
and  beam  floor  were  tested  to  failure  under  static 
uniform  load  on  the  interior  panel.  Tests  were 
made  at  the  University  of  Bristol.  Panels  were  12in. 
square  and  varied  in  depth  from  about  3/8in.  to 
about  ll/l6in.  Twenty  slabs  were  tested,  varying 
in  ages  from  9  to  14  days.  High  abutment  stresses 
were  manifested  by  diagonal  cracks  across  a  corner 
of  the  surround,  consisting  of  the  eight  exterior 
panels.  The  ultimate  loads  of  the  interior  panels 
of  slabs  with  continuous  tee  reinforcement  showed 
good  agreement  with  theory  which  includes  compressive 
membrane  strains  in  the  panel,  and  lateral  edge 
displacements  once  the  concrete  had  cracked,  but 
did  not  reach  the  results  of  theory  which  ignored 
these  strains  and  displacements.  Results  emphasize 
the  importance  of  abutment  stiffness  in  the  mobili¬ 
zation  of  arching  forces. 

4.  Christiansen, K.  P. ,  The  Effect  of  Membrane  Stresses 
on  the  Ultimate  Strength  of  Interior  Panel  in  a 
Reinforced  Concrete  Slab,  The  Structural  Engineer 
(London)  Vol.  41,  No.  8,  August  196 3". 

Christiansen,  a  Civil  Engineer  in  Copenhagen,  per¬ 
formed  the  reported  work  and  wrote  the  paper  while 
a  lecturer  at  the  University  of  Wltwatersrand, 
Johannesburg,  South  Africa.  The  author  presents 
his  theory  for  evaluation  of  arching  action  in  slabs 
and  reports  results  of  tests  of  4  identical  pairs 
of  simply  supported  beams,  with  one  of  each  pair 
axially  restrained  by  a  welded  steel  frame  and  the 
other,  unrestrained;  all  beams  were  6  Inches  wide 
and  spans  and  depths  were  60  x  3*5  Inches,  60  x  3 
inches  and  '(?  x  3  inches.  Load  carried  by  arching 
In  the  axially  restrained  beams  was  30  -  35*  of 
the  total  In  three  cases  with  cube  strengths  from 
4090  to  4975  ns l,  and  705  of  the  total  in  a  60  x  3*5 
Inch  beam  with  cube  strength  of  5660  psi.  Author 
concludes  that  although  loads  supported  by  arching 
can  be  predicted  with  some  accuracy,  more  research 
is  required  to  make  arching  action  a  reliable  and 
useful  part  of  design. 
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5.  Self,  W.  M. ,  Ultimate  Strength  of  Reinforced  Concrete 
Flat  Slabs,  Journal,  Structural  Division,  ASCE 

Vol .  90,  No.  St 4,  August  1964,  P.205.  Three 
reinforced  concrete  slab  models  were  tested.  Slabs 
were  10  feet  square  by  2  inches  thick,  and  supported 
by  four  interior  columns  producing  an  interior 
panel  5?s  feet  square.  Overhang  balanced  moments 
over  columns,  thus  eliminating  influence  of  column 
stiffness . 

Yield  line  pattern  conformed  to  predictions  in 
accordance  with  Johansen  theory.  Ultimate  strength 
of  slab  No.  1,  reinforced  with  cold-rolled  steel, 
was  9%  above  prediction;  slab  No.  2,  reinforced  with 
hot-rolled  steel,  reached  an  ultimate  strength  17% 
above  prediction.  (Slab  No.  3  failed  in  shear  before 
flexural  yield).  The  variation  is  indicative  of 
strain  hardening  of  steel.  The  slabs  lacked  sur¬ 
rounding  confinement  necessary  to  mobilize  compres¬ 
sive  membrane  forces. 

6.  Hodge,  P.  G.,  Jr.  and  Perrone,  N.,  Yield  Load  of 
Slabs  with  Reinforced  Cutouts,  Journal  of  Applied 
Mechanics,  Vol.  24,  March  1957,  PP. 85-92. 

Presents  a  method  of  computing  the  upper  and  lower 
bounds  on  the  capacity  of  a  plane  slab  with  a 
reinforced  cutout  to  carry  uniaxial  tensile  loads 
which  are  either  uniform  or  applied  by  means  of  a 
perfectly  rigid  clamp.  Results  are  applied  to 
several  examples  and  found  to  agree  quite  well 
with  experimental  values. 

7.  Gamble,  W.  L. ,  Sozen,  M.  A.,  and,  Siess,  C.  P. , 
Measured  and  Theoretical  Bending  Moments  in  Rein¬ 
forced  Concrete  Floor  Slabs,  University  of  Illinois 
Civil  Engineering  Studies,  Structural  Research 
Series  No.  246,  June  1962. 

Test  results  of  two-way  slabs  reported  in  the  lit¬ 
erature  and  results  of  five  tests  of  differing  1/4- 
scale  slab  models,  performed  at  the  University  of 
Illinois,  were  reviewed  with  the  objective  of 
developing  a  slab  design  procedure. 

Slab  strengths  were  calculated  according  to  yield 
line  theory,  Moe's  formula  for  shear  strength, 
Elstner  and  liognestad  punching  shear  equations  or 
the  equation  from  ACI  committee  326. 
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In  nearly  all  of  the  tests  reported  in  the  lit¬ 
erature,  the  specimens  failed  at  loadings  higher 
than  cal  alated  according  to  yield  line  theory. 

Some  reasons  advanced  for  the  low  calculated  values 
are : 

1.  Incorrect  determination  of  steel  strength. 

2.  Strain  hardening  of  reinforcement. 

3.  Errors  in  establishing  location  of  yield 
line . 

4.  Arching  of  compressive  membrane  forces. 

5.  Effects  of  deformation  of  the  structure. 

In  most  of  the  five  University  of  Illinois  tests, 
the  flexural  mode  of  failure  was  complicated  by 
other  distress  (shear  failure,  beam-column  distress, 
etc.).  The  test  results  indicated  that  if  flexural 
failure  had  been  allowed  to  occur,  the  failure 
load  would  have  exceeded  the  results  obtained  by 
yield  line  analysis. 

Hillerborg,  A.,  A  Plastic  Theory  for  the  Design  of 
Reinforced  Concrete  Slabs ,  Preliminary  Publication, 
International  Association  for  Bridge  and  Structural 
Engineering,  Sixth  Congress,  Stockholm,  I960. 

Conclusion  is  reached  that  although  the  yield  line 
theory  will  reveal  the  magnitude  of  a  load  great 
enough  to  cause  failure  through  the  formation  of 
plastic  hinges,  the  load  is  theoretically  unsafe 
since  other  yield  lines  may  form  at  a  lower  level. 

The  equilibrium  theory  states  that  if  a  distribu¬ 
tion  of  moments  can  be  found  which  satisfies  the 
equilibrium  equation  and  the  boundary  conditions, 
under  the  action  of  a  given  load,  and  these  moments 
do  not;  exceed  the  yield  moments  at  any  section  of 
the  slab,  the  slab  is  capable  of  carrying  that  load. 
This  is  a  lower  bound  and  the  exact  ultimate  load  is 
somewhere  between  the  results  of  the  yield  line  and 
equilibrium  theories. 

Hillerborg,  A.,  Theory  of  Equilibrium  for  Rein¬ 
force  l  Concrete  Slabs,  Department  of  Scientific 
and  Industrial  Research,  Building  Research  Station, 
Library  Communication  No.  1082,  Great  Britain. 
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If,  for  a  certain  external  load,  q(x»y)»  a  moment 
distribution  can  be  found  which  satisfies  the 
equation , 


52>>1x  +  *2mi<  *2"xtf 

3X2  3</2  3x3tf 


and  the  edge  conditions,  and  if  the  slab  can  take 
up  these  moments  at  each  point,  q(x,y)  is  a  lower 
limiting  value  for  the  bearing  capacity  of  the  slab. 
The  moments  are  not  coupled  and  therefore,  two 
may  be  selected  arbitrarily  and  the  remaining  one 
solved. 

A  practical  solution  is  advanced,  involving  the  di¬ 
vision  of  the  slab  into  strips,  in  which  the  load 
in  one  strip  is  carried  in  one  direction  only. 

10.  Crawford,  Robert  E. ,  Limited  Design  of  Reinforced 
Concrete  Slabs ,  Journal  of  the  Engineering  Mechanics 
Division,  American  Society  of  Civil  Engineers, 

Vol.  90,  No.  EMS,  October,  1964. 

The  yield  line  theory  is  considered  in  terms  of 
limit  analysis  and  Is  shown  to  give  an  upper  bound 
on  the  collapse  load.  Similarly,  the  equilibrim 
theory  is  reviewed  and  shown  to  give  a  lower  bound. 

11.  Kemp,  K.  0.,  A  Lower  Bound  Solution  to  the  Collapse 
of  an  Orthotropioally  Reinforced  Slab  on  Simple 
Supports ,  Magazine  of  Concrete  Research  (London) 
July,  1962. 

A  lower  bound  solution  to  the  collapse  of  a  simply 
supported,  rectangular  slab,  orthotropically  rein¬ 
forced  and  carrying  a  uniformly  distributed  load  is 
developed.  The  lower  bound  collapse  loads,  cal¬ 
culated  for  a  range  of  coefficients  of  orthotropy 
and  ratio  of  width  to  length  of  slab,  agree  closely 
with  the  upper  bound  values  derived  from  the  yield 
line  theory.  The  yield  lines,  however,  are  quite 
different;  they  are  found  to  be  an  infinite  set  of 
curves  in  the  lower  bound  solution.  With  the  solu¬ 
tion,  the  negative  reinforcement  required  in  the 
corners  of  the  slab  and  the  loads  transmitted  to 
the  supporting  edges  were  determined. 
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12.  Brotchie,  John  F.  ,  Jacobson,  Amnon  and  Okubo,  Sadaji, 
Effect  of  Membrane  Action  on  Slab  Behavior ,  Report , 
Department  of  Civil  Engineering,  M.I.T.,  U.S.  Naval 
Civil  Engineering  Laboratory,  Port  Hueneme,  Calif¬ 
ornia,  August,  1965. 

Forty-five  square  slabs  of  plain  and  reinforced 
concrete  were  tested  under  uniform  loading.  All 
spans  were  15"  x  15".  Thicknesses  were  0.75", 

1.5"  and  3*0",  resulting  in  span/depths  ratios  of 
20,  10  and  5. 

Lower  reinforcement  only  was  used,  distributed 
uniformly  and  equally  in  each  direction.  Reinforce¬ 
ment  ratios  were  0% ,  0.5%,  b% ,  2%  and  3%. 

Five  sets  of  boundary  conditions  were  used: 

1.  Restrained  at  edges  against  axial  elonga¬ 
tion  only,  at  approximate  level  of  rein¬ 
forcement,  by  24  cells. 

2.  Same  as  1,  with  added  resistance  to  internal 
shear  in  the  slab  at  the  support. 

3.  Same  as  1,  but  with  level  of  restraining 
force  raised  to  the  middle  of  the  edge 
surface  of  the  slab. 

4.  Slab  clamped  at  the  supports  by  top  plate 
and  base  plate  and  restrained  at  the 
edges  by  epoxy  resin  fill. 

5.  Simply  supported  on  0.75"  diameter  roller 
bearings . 

The  steel  restraining  frame  was  designed  for  essen¬ 
tially  complete  lateral  and  vertical  rigidity.  Twelve 
of  the  24  restraining  cells,  uniformly  distributed 
along  two  adjacent  edges,  were  wired  with  electrical 
resistance  strain  gauges  for  measuring  the  restraining 
force . 

The  effect  of  arching  was  found  to  be  significant  in 
slabs  with  span/depth  ratios  of  20;  equivalent  in  load 
capacity  to  approximately  2%  of  conventional  rein¬ 
forcement.  The  effect  of  arching  is  greater  for  thick 
slabs  and,  for  span/depth  ratios  of  5,  is  equivalent 
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to  over  31  of  conventional  reinforcement.  External 
restraint,  mobilizing  arching  action,  increases  stiff¬ 
ness  and  reduces  cracking.  When  maximum  arching  - 
supported  load  is  reached,  however,  load  capacity  drops 
while  deflection  increases,  allowing  sudden  failure  if 
the  slab  is  lightly  reinforced  and  the  applied  loading 
is  sustained. 

The  magnitude  and  distribution  of  the  restraining  force 
indicates  that  essentially  the  full  capacity  of  the 
concrete  in  compression  is  utilized,  when  the  slab  is 
restrained. 

For  thin  slabs,  essentially  full  edge  restraint  is 
necessary  to  mobilize  the  full  increase  in  load  capa¬ 
city.  For  thicker  slabs,  additional  edge  displacement 
may  be  tolerated,  without  significantly  reducing  the 
load  capacity. 

The  effect  of  direct  normal  stress  due  to  the  uniform 
load,  in  combination  with  two-way  inplane  stresses, 
is  to  produce  a  triaxially  compressed  zone.  In  the 
case  of  very  thick  slabs,  that  is,  span/depth  ratios 
of  5j  the  effect  is  apparently  significant;  the  plastic 
deformability  of  the  concrete  in  the  compressive  zone 
is  notably  increased. 

Tensile  membrane  action  is  significant  only  at  deflec¬ 
tions  which  are  too  large  for  use  in  resisting  service 
loads.  The  tensile  membrane  action  may  have  some 
appeal,  however,  to  provide  protection  against  col¬ 
lapse  in  the  case  of  a  single  catastrophic  loading, 
particularly  in  view  of  the  energy  absorption  which 
precedes  in  this  range. 

Simple  theoretical  expressions  are  presented  which 
were  found  to  predict,  with  some  exceptions,  the  ranges 
of  behavior  described. 
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Appendix  A 


COMPUTER  SUBROUTINE  FOR  CALCULATING  LOADS 
OH  AN  UNDERGROUND  BOX 

The  formulation  of  the  load  subroutine  contained  in 
the  computer  nroqram  is  described  in  this  aopendix. 
Pressures  felt  by  a  shallow  buried  box  are  calculated 
usinq  the  enninecrinp;  approach  of  Reference  2  (Chapter  iJ), 
Soil-structure  interaction  has  not  been  considered,  ner 
se.  The  loading  formulation  assumes  that  the  overpressure 
of  interest  are  less  than  100  psi.  Also  megaton  yield 
weapons  are  assumed. 

A-l  Airblast  Characteristics 

This  section  is  a  brief  statement  of  the  airblast 
characteristics,  as  recorded  by  Brode  in  /  Review  of 
Nuclear  Explosion  Phenomena  Pertinent  to  Protective 
Construction** ,  This  reference  should  be  consulted  for 
a  complete  review  of  the  physical  phenomena  pertinent  to 
protective  construction. 

Figure  A-l  shows  the  shock  arrival  time  ts  and  the 
shock  radius  Rs  for  the  overpressures  resulting 
from  the  detonation  of  a  20  H?  bomb.  These  values 
depend  upon  the  energy  of  the  explosion  and  may  be 
scaled  by  the  cube  root  of  the  yield  ratio  in  mega¬ 
tons  for  yields  other  than  20  MT.  This  figure  also 
illustrates  the  shock  velocity  as  a  function  of  the 
peak  overpressure  which  is  independent  of  weapon 
yield. 

The  time  history  of  the  overpressure  is  described 
quite  well  at  all  pressure  levels  under  100  psi  by 
the  sum  of  two  decreasing  exponential  functions  of 
time : 

tP  »  6P  (Ae‘aT  ♦  bc'eT)(I-x)  (A-l) 

8 

where  t  is  the  time  after  shock  arrival  measured 
in  units  of  positive  phase  duration.  To  force  this 
curve  to  to  zero  overpressure  at  the  end  of  the 
positive  phase,  a  linear  factor  has  bet..i  included 
that  becomes  zero  at  a  tine  equal  to  the  duration 
of  the  positive  phase  (t *1  when  “POt 

8  p 
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Figure  A-l  SHOCK  RADIUS  (Ks),  ARRIVAL  TIME  U3) 
AND  SliOCK  STEED  ((/)  VERSUS  PEAK  OVERPRESSURE 
UP5)  FOR  A  2 OMT  SURFACE  BURST.  (Ref.  15) 


Figure  A-2  gives  the  values  of  all  shock  parameters 
and  coefficients  necessary  to  obtain  the  pressure¬ 
time  curve  for  a  given  peak  overpressure. 

Curves  showing  the  pressure  time  relations  based  on 
this  analytical  expression  are  given  in  Figure  A-3. 

A-2  General  Characteristics  of  the 
Soil  Pressure  Wave 


The  alteration 
pressure  wave  as  it 
trated  in  Figure  A-Jl 
very  rapid  rise,  an 
time,  Pp,  similar  to 
At  greater  deaths  a 
somewhat  lengthened 
lengthened  rise  time 
The  impulse  --  total 
constant. 


in  free-fjeld  form  of  the  vertical 
progresses  into  the  soil  is  illus- 
The  surface  soil  pressure  has  a 
exponential  decay,  and  a  duration 
that  of  the  airblast  pressure  wave, 
reduced  soil  pressure  oeak  and  a 
duration  are  Dresent.  Also,  a 
to  peak  soil  pressure  takes  place, 
area  under  the  curve  —  is  nearly 


The  attenuation  of  stress  with  depth  is  generally 
attributed  to  two  causes,  energy  absorotior  by  the  soil 
and  by  three  dimensional  dispersion  of  the  energy.  In¬ 
formation  regarding  energy  absorption'  is  not  available 
In  a  form  useful  for  calculating  loads.  Three-dimension¬ 
al  or  spatial  attenuation  can  be  aooroximated  by  the 
following  expression  proposed  by  Newmark2: 
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where 


b?_ 


VPt 


u 


peak  overpressuee  at  the 
surface,  psl 

maximum  vertical  soil  stress 
at  depth  n  (ft),  psi 

geometric  attenuation  factor 


W  ■  yield  in  megatons 


Figure  A-3  FORM  FACTORS  FOR 
OVERPRESSURES  FROM  NUCLEAR  EXPLOSIONS 

(Ref.  15) 
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The  above  approximation  does  not  permit  soil  stress 
higher  than  the  peak  overpressure,  a  condition  that  can 
truly  apply  only  when  the  airblast  velocity  is  low  com¬ 
pared  to  the  seismic  velocity, lb 

At  the  100  psi  range  of  a  20  MT  surface  burst,  the' 
geometric  attenuation  factor  (Equation  A-3)  for  100  feet 
of  depth  is  0,86,  For  lower  overpressures,  higher  yields 
or  less  depth,  the  attenuation  factor  approaches  1,0, 

Referring  again  to  Figure  A-*J ,  it  is  obvious  that 
the  time  required  for  the  vertical  soil  pressure  to 
reach  its  peak  value,  can  never  be  less  than  the  air- 
blast  rise  time.  Figure  2-2.3  of  Reference  3  gives  an 
indication  of  the  magnitude  of  this  airblast  rise  time, 
as  it  varies  with  overpressure.  However,  the  value  of 
rise  time  is  not  well  defined. 

Reference  2  expresses  the  rise  time  in  soil  as; 
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( A-5) 


where 


r 


* rab 


airblast  rise  time 


y  ~  depth  of  interest 


seismic  velocity  consistent  with 
the  loading  modulus  of  the  soil 


c  .  =  seismic  velocity  consistent  with 
the  initial  tangent  modulus  of 
the  soil. 


When  the  values  of  ^p  and  c-i  are  not  well  defined, 
an  approximate  evaluation  of  the  rise  time  must  be  made. 

In  view  of  the  many  uncertainties  regarding  airblast 
rise  time  and  soil  properties,  Reference  2  suggests  an 
approximation;  that  the  rise  time  be  taken  as  one-half 
the  time  required  for  the  peak  stress  to  reach  the  depth 
of  interest, 

t  -  lJ-  (A-6) 

This  expression  (Equation  A-6)  has  been  used  within  \ 
the  computer  program. 
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Horizontal  stress  in  the  soil,  P^,  is  taken  as  some 
constant  times  the  vertical  soil  stress, 

Ph  -  *  Pv  (A-7) 

The  value  of  the  constant  K,  depends  upon  the  soil  prop¬ 
erties,  the  degree  of  saturation,  the  stress  level,  and 
the  condition  of  lateral  restraint  at  the  depth  of  in¬ 
terest. 


A-3  Loads  on  the  Box 

Each  finite  element  of  the  box-structure’s  exterior 
roof  and  walls  is  assigned  a  load-time  history.  Figure 
A- 5  shows  a  roof  element,  a  side  wall  element,  and  a 
front  wall  element.  This  figure  gives  each  dimension 
used  in  calculating  the  idealized  pressures  on  the  box. 
Table  A-l  lists  the  formulas  used  in  calculating  the 
form  of  the  pressure  on  each  finite  element,  expressed 
as  an  average  uniform  pressure  over  the  element, 

A-3.1  Roof  Pressure 


The  vertical  Dressure  on  a  roof  element  initiates 
with  the  arrival  of  the  shock  front  at  the  leading  edge 
and  builds  up  to  a  maximum  as  the  shock  front  reaches 
the  trailing  edge.  The  average  pressure  is  assumed  to 
increase  linearly  from  zero  to  its  maximum  value  during 
the  time  required  for  the  shock  front  to  travel  across 
the  element  (see  Table  A-l), 

Since  the  box  roof  is  divided  into  a  number  of 
finite  elements  for  analysis  purposes,  the  average  pres¬ 
sure  is  considered  element-by-element  in  the  direction 
of  the  shock  travel. 

Arching  of  the  soil  above  the  roof  is  not  considered. 
Also,  spatial  attenuation  of  the  peak  pressure  with  depth 
of  cover  is  not  considered, 

A-3. 2  Wall  Pressure 


The  horizontal  pressures  in  the  soil  are  calculated 
empirically  by  multiplying  the  vertical  pressures  by  the 
coefficient  of  earth  pressure  at  rest,  K0,  The  analyst 
choses  the  value  of  K 0  believed  appropriate  for  the 
added  increment  of  airblast  induced  horizontal  stress. 
Suggested  values  are  given  in  Table  A-2, 
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TABLE  A-l 

AVERAGE  PRESSURE  ON  AN  ELEMENT 
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No  spatial  attenuation  of  the  peak  horizontal  soil 
pressure  has  been  included  in  the  programed  calcula¬ 
tions.  Very  little  spatial  attenuation  is  predicted  by 
Equation  A-3  for  depths  less  than  100-ft.,  located 
within  the  100  psi  range  of  megaton  weapons.  The  co¬ 
efficients  of  lateral  earth  pressure  at  rest  (Table 
A-2)  are  too  crude  to  be  modified  further  by  a  depth 
attenuation  factor. 

The  horizontal  pressure  on  a  wall  element  is  as¬ 
sumed  to  initiate  with  the  arrival  of  a  vertical  com¬ 
pression  wave  at  the  top  edge,  which  is  induced  by  the 
shock  front  as  it  arrives  directly  overhead.  The 
average  pressure  is  assumed  to  increase  linear? y  from 
zero  at  initial  arrival  to  a  maximum  value  of  ^oApa 
during  the  time  required  for  the  soil  pressure  pulse  to 
travel  across  the  element  plus  the  rise  time  calculated 
from  Equation  A-6,  (See  Table  A-l.) 

Since  the  box  is  divided  into  a  number  of  finite 
elements,  the  average  pressure  is  considered  element- 
by-element  with  increasing  depth  and  in  the  direction 
of  shock  travel. 

A-3. 3  Boundary  Conditions 

As  the  shock  front  passes  over  the  box,  the  sur¬ 
rounding  ground  surface  displaces  vertically  (Figure 
A-6).  Overpressure  loads  are  transmitted  vertically 
through  the  box  to  the  foundation  material  immediately 
below.  A  distributed  elastic  spring  is  used  to  repre¬ 
sent  the  penetration  compliance  of  the  foundation.  Thus, 
the  vertical  force  system  is  one  of  overpressure  forces 
driving  the  box  into  the  foundation  below.  Frictional 
forces  on  the  walls  are  ignored.  The  value  of  the  found¬ 
ation’s  penetration  stiffness,  K  ,  is  chosen  by  the  ana¬ 
lyst. 

Horizontal  pressures  act  concurrently  with  the  verti¬ 
cal  pressures.  A  horizontal  soil  spring  value  of  stiff¬ 
ness,  Ku,  is  chosen  by  the  analyst.  The  program  converts 
the  average  element  pressures  into  displacements  at  the 
far  end  of  this  horizontal  soil  spring  by  dividing  the 
pressures  by  the  spring  stiffness  (see  Figure  A-7).  The 
resulting  time-dependent  displacement  is  treated  as  the 
horizontal  forcing  function. 
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Figure  A-6  BOUNDARY  CONDITIONS, 
VERTICAL  RESPONSE. 


Figure  A- 7  BOUNDARY  CONDITIONS, 
HORIZONTAL  RESPONSE 
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TABLE  A- 2 


Ratio  of  Horizontal  to  Vertical  Soil  Pressures 

(Reference  2) 


K0 ,  For  Stresses  Up  to  1,000  psi 


Soil  Description 

Dynamic 

Static 

Undrained 

Undrained 

Drained 

Cohesionless  Soils, 

Damp  or  Dry 

1/4 

1/3-dense 

1/2-loose 

1/3-dense 

1/2-loose 

Unsaturated  Cohesive 
Soils  of  Very  Stiff  to 
Hard  Consistency 

1/3 

1/2 

1/2 

Unsaturated  Cohesive 
Soils  of  Medium  to 

Stiff  Consistency 

1/2 

1/2 

1/2 

Unsaturated  Cohesive 
Soils  of  Soft 
Consistency 

-3/4 

1/2  to  3/4 

1/2  to  3/4 

Saturated  Soils  of 

Very  Soft  to  Hard  Con¬ 
sistency  and  Cohesion¬ 
less  Soils 

1 

1 

1/2-stiff 

3/4-soft 

Saturated  Soils  of  Hard 

Consistency, 

qu  •  4  tsf  to  20  tsf . 

3/4  to  1 

1 

1/2 

Saturated  Soils  of  Very 
Hard  Consistency. 
qu  ■  20  tsf. 

3/4 

1 

1/2 

Rock  Obtain  from  tests  on  rock 

-  cores  and  correlate  with 
seismic  data. 
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A-3.*J  Outrunning  Ground  Actions 


Many  structures,  which  are  designed  for  overpres¬ 
sures  less  than  100  psi,  will  be  in  the  region  of  out¬ 
running  ground  motion.  The  phasing  between  outrunning 
ground  notions  and  airblast  arrival  can  be  calculated 
by  the  procedure  described  in  Reference  8  and  need  not 
be  repeated  here.  The  phasing  time  is  part  of  the  in¬ 
put  data  furnished  by  the  analyst. 

This  section  is  a  statement  of  ground  motion  pre¬ 
diction  as  given  by  Sauer.  The  equations  used  in  the 

computer  program  to  describe  the  amplitude  of  outrun¬ 
ning  ground  notion  are  presented.  In  addition,  the 
particle  velocity  waveform  is  presented  for  both  verti¬ 
cal  and  horizontal  motion. 

.  The  vertical  velocity  waveform  (Figure  A-8)  is 
completely  identified  by: 

1^2  ,  the  amplitude  of  the  third 
velocity  peak 


t.2 ,  the  period  of  the  first 
three  half  waves 


These  variables  are  calculated  from  the  following 
relationships : 


4x10n 
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^(wsec)  **  100  +  AR/4 


- (ft/sec) 


(A-8) 
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where, 

aP8  ■  Peak  overpressure  (psi) 

»  Range  from  ground  zero  (ft.) 

&p  «  Distance  from  site  to  point  of  initial  out¬ 
running  (ft.) 

a  • In  situ  specific  gravity 
«»  Adjusted  seismic  velocity. 

V *  »  for  soils  and  incompetent  rock; 

*  VL  for  competent  rock, 

m  Seismic  velocity  (ft/sec) 


W  ■  Weapon  yield  (negatons). 
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Figure  A- 8  NORMALIZED  VERTICAL  WAVEFORM 
OUTRUNNING  GROUND  MOTION  (Ref. 8) 


Figure  A-9  RECOMMENDED  PROCEDURE  FOR 
HORIZONTAL  OUTRUNNING  WAVEFORM  (Ref.  16) 
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The  displacement  waveform  of  Figure  A-8,  Pv(-t),  is 
stored  in  the  program.  This  free-field  displacement  is 
applied  through  the  foundation  springs  (Figure  A-6) 
when  the  analyst  enters  the  above  data  on  the  input 
forms.  For  computational  purposes  the  equation 
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*max 
2.5  tz 


*wax  | 
2.5  .J 


(A-10) 


is  adopted  to  represent  the  horizontal  waveform  shown  in 
Figure  A-9. 


For  t  ■  0  and  t  *  2.5  *2,  l/2(0)  ■  0. 


2.5  f2  is  the  duration  of  the  wave,  K  fixes  the  ampli¬ 
tude,  and  m  and  p  are  constants  responsible  for  the  shape 
of  the  curve.  Quantity  m  is  related  to  p  by 


m  «  p 


wax 

2.5  t. 


wax 

2.5  t , 
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The  program  Integrates  Equation  A-10.  giving  the  free 
field  horizontal  soil  displacements,  P.  C I .  The  wall  load¬ 
ing  used  in  the  computer  program  assumes  that  the  free- 
field  displacements,  P,(X|,  are  applied  through  the  effec¬ 
tive  soil  springs  (Figure  A-7). 
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»  ABSTRACT 


this  report  describes  the  development  of  two  digital  computer 
programs.  Both  programs  analyse  the  response  of  a  multi-cell 
rectangular  box-structure  to  external  dynamic  loads.  The  programs 
use  the  principle  of  finite  elements  to  represent  the  box-structure 
and  solve  the  equations  of  equilibrium  by  mesne  of  e  Gaussian 
elimination  procedure.  The  first  of  the  two  programs  considers 
the  dynamic  response  of  the  box-structure  with  Its  material 
properties  remaining  linear.  The  output  consists  of  displacements, 
moments  and  stresses.  Also,  response  spectra  may  be  generated  at 
equipment  attachment  points  within  the  structure.  The  second 
program  calculates  the  dynamic  response  with  the  added  provision 
that  when  moments  In  any  plate  exceed  a  predetermined  yield  value, 
that  plate  will  not  accept  any  increase  in  loads  normal  to  Its 
plane.  However,  It  Is  assumed  that  the  plate  would  still  be  able 
to  carry  the  Inplane  diaphragm  forces. 

A  synoptic  review  of  selected  literature  pertaining  to  the  failure 
modes  and  failure  loadings  of  two-way  slabs  Is  presented  at  the 
end  of  the  report. 
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